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ABSTRACT 



Nuclear physics can be applied in various ways to the study of neutron stars. This thesis reports 
on one such application, where the relativistic mean-field approximation has been employed 
to calculate the equations of state of matter in the neutron star interior. In particular the 
equations of state of nuclear and neutron star matter of the NL3, PK1 and FSUGold parameter 
sets were derived. A survey of available literature on neutron stars is presented and we use the 
derived equations of state to reproduce the properties of saturated nuclear matter as well as the 
mass-radius relationship of a static, spherical symmetric neutron star. Results are compared 
to published values of the properties of saturated nuclear matter and to available observational 
data of the mass-radius relationship of neutron stars. 
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OPSOMMING 



Kernfisika kan op vele maniere aangewend word binne die studie van neutronsterre. Een so 'n 
toepassing is die gebruik van die relatiwistiese gemiddelde-veld benadering om die toestandsverge- 
lyking van neutronstermaterie af te lei. Die afleiding van die toestandsvergelyking van kern- en 
neutronstermaterie vir onderskeidelik die NL3, PK1 en FSUGold parameter stelle vorm die ba- 
sis van die tesis. 'n Oorsig oor die beskikbare literatuur aangaande neutronsterre word gebied 
en van die gepubliseerde resultate word herbereken. In die besonder word die genoemde toe- 
standsvergelykings gebruik om die eienskappe van versadigde kernmaterie te bepaal, asook om 
die massa-radius verhoudings van statiese, sferies-simmetriese neutronsterre mee te bereken. Die 
resultate word vergelyk met gepubliseerde waardes vir die eienskappe van versadigde kernmaterie 
en die waargenome massa-radius verhoudings van neutron sterre. 
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CHAPTER 1 
Introduction 



Judging from current experimental activities, one of the pertinent questions that nuclear physicists 
are trying to answer is: What was the state of matter, at a time when the universe was very 
hot and dense? High-energy collisions of heavy- ions, such as the ones performed in the PHENIX 
experiment at the Relativistic Heavy-Ion Collider (RHIC) and planned at CERN in the ALICE 
experiment, aim to study such hot, dense nuclear matter [U[2]- More specifically the experiments 
at RHIC aim to study the state of matter, called quark-gluon plasma, which is believed to be 
the state of matter shortly after the Big Bang [TJ. 

Another question, that is currently investigated in earnest by nuclear physicists, is if the structure 
of normal nuclear matter (protons and neutrons) breaks down under conditions of low tempera- 
ture and extreme pressure, i.e. whether the ground state of nuclear matter consists out of matter 
other than nuclear matter. To answer this question, nuclear physicists must look further afield 
than laboratory experiments. In laboratory experiments, where either two beams of high energy 
particles are collided or high energy particles are shot onto a target, a highly energetic, but short- 
lived state of matter is created, which is rather far removed from the ground state of matter [3]. 
A physical system is in its ground state when the system is in its lowest energy configuration. 
The ground state is usually achieved after the system had some time to equilibrate. 56 Fe is the 
nucleus with the greatest (negative) binding energy and therefore has the lowest ground state of 
all the nuclei. Questions have been raised whether 56 Fe is the ultimate ground state of matter 
since W. Baade and E. Zwicky proposed in 1934 that some supernovae are driven by the binding 
energy of a neutron star [3], i.e. that a massive star explodes due to the release of energy in 
forming a neutron star. On average neutron stars have a mass of about 1.4 times the mass of 
our sun [5]. The gravitational binding energy of a neutron star is about 10% of its mass, while 
the nuclear binding energy of 56 Fe is about 9 MeV/nucleon, which is about 1% of the mass of 
a 56 Fe nucleus [4]. Thus the gravitational binding energy of a neutron star consisting of 56 Fe 
would be about ten times larger than the nuclear binding energy. Since the star is in a stable 
long-lived state, one can assume that the matter is in a state other than that of 56 Fe. As such 
nature has provided us with a laboratory to study cold, dense matter, which we currently cannot 
synthesise in the laboratory. (By cold is meant that the thermal energy of the particles is very 
small compared to the fermi energy and thus thermal excitations are assumed not to take place 

IS!-) 
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Neutron stars are believed to be formed in the core-collapse supernovae of massive stars [HE]- A 
normal star is stabilised against gravitational collapse by the thermal pressure due to the energy 
release of nuclear fusion processes in the star. Fusion in the core of the star would proceed the 
fastest due to the higher pressure at the centre of the star. Once these fusion processes have 
reached the formation of 56 Fe, the last exothermic fusion phases, it will shut down and the core 
will start to cool and contract under its own gravity [4]. As the core contracts densities where 
electrons become relativistic will quite easily be reached. The energy of the core can thus be 
reduced by the capture of electrons by the protons (inverse beta-decay), thus making the core 
more neutron-rich [7]. The crushing effect of gravity will be halted by the short-ranged repulsion 
of the strong nuclear force. 

Matter falling onto the core, due to the low pressure created by the contracting core, would 
rebound off the stiffened core, creating a shockwave that travels outward from the centre of the 
star, which stalls after some hundred kilometres, due to energy losses as the wave travels through 
the interior of the star. This stalled shockwave creates an accreting front as matter from the 
outer parts of the star collapses towards the centre of the star. Through poorly understood 
mechanisms the binding energy of the neutron star (the core of the collapsing star) gets trans- 
ferred to the accreting shock front, which causes the rest of the star to explode in a supernova. [1] 

In 1967 Jocelyn Bell observed a pulsating radio-source in outer space that had characteristics 
unlike any other radio-source [9j. Initially the origin of these radio-pulses was unclear (little 
green men were not ruled out), but eventually the source was explained to be a rapidly rotating 
neutron star, today known as a pulsar |10t lllj. Today more than 1100 pulsars, the commonly 
observed state of neutron stars, are known [4]. 

As with the question regarding the ultimate ground state of matter, the constituents of the 
interior of neutron stars are also not known. As the average densities of neutron stars are com- 
parable to that of nuclei, it is assumed that neutron stars consist (at least in some part) of 
baryonic matter (such as protons and neutrons) and therefore can be viewed as giant nuclei, 
but with a mass number of 10 57 [7\\ The main difference between nuclei and neutron stars is 
that while nuclei are bound by the nuclear strong force, neutron stars are bound by gravity. 
The fact that gravity is attractive on all scales implies that the neutron star must have some 
form of internal pressure to counteract the effect of gravity, otherwise no stable neutron stars 
would exists, only black holes (objects that have collapsed under their own gravity). Baade and 
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Figure 1.1: A representation of different ideas regarding the composition of the 
interior of neutron stars, from Ref. |8j 

Zwicky envisioned that a neutron star is supported against gravitational collapse by the nucleon 
degeneracy pressure: the pressure that is due to Pauli's Exclusion Principle that states that no 
two identical particles can occupy the same energy state. Thus the pressure is given by particles 
that all want to occupy the lowest energy state, but these states are filled from the bottom and 
once they are full particles have to occupy higher lying states [4]. Due to the central role of 
gravity, neutron stars can only be adequately described using the general theory of relativity [4j. 
In Chapter [2] necessary concepts in relativity will be introduced to formulate such a description, 
which will be derived in Chapter [3j 

But that would not be the full story. To be able to have such an adequate, relativistic description 
of neutron stars information regarding the relation between the pressure and the energy density 
(equation of state) of matter in the interior of the neutron star is needed. If we assume that the 
neutron star consists of baryonic matter, models of nuclear matter can be used to provide the 
equation of state of the neutron star interior [U EJ [7] . Fig. 11.11 shows a schematic presentation 
of various assumptions about the interior of neutron stars (and as such for the ground state of 
matter). 

These different models for the interior of the neutron star can be tested by comparing calculated 
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properties, specifically the mass-radius relationship fo a neutron star, to the observed values. By 
this process suitable models (and equations of state derived from these models) can be identi- 
fied for the description of cold, dense matter. In general these different models of nuclei and 
nuclear matter are complex and difficult to solve exactly and therefore these models have to be 
approximated. The focus of this work will be to derive the equation of state of neutron star mat- 
ter from a specific model of nuclei and nuclear matter, called quantum hadrodynamics (QHD), 
through the application of relativistic mean-field theory. This will be done by making a survey 
of some of the available literature on this subject as well as to try to reproduce some of the 
published values of certain properties of nuclear matter and neutron stars. The aim of this work 
is to consolidate some of the current knowledge in the theoretical study of neutron stars. Since 
these theoretical studies can only be validated through good agreement with observed properties 
of neutron stars, close co-operation between theorists and astronomers are crucial in the study of 
neutron stars. The converse is also true: to explain observational results theoretical modelling is 
needed to predict certain observed properties. Only through this interplay between theoretical 
and observational science can the understanding of our universe be advanced. This work there- 
fore also aims to include some of the observational results to constrain the theoretical calculations. 

It is foreseen that much more will be learned about neutron stars with the commission of the 
Square Kilometre Array (SKA). The SKA will be the world's biggest radio-telescope, with an 
effective collecting area of one square kilometre, which will make the SKA fifty times more 
powerful than the current most powerful radio-telescope |12l 13J. To explain the observations 
made by the SKA knowledge of the theoretical modelling of neutron stars would be crucial. 

1.1 South African perspective 

This work is applicable within the South African context since South Africa is currently 
developing the meerKAT radio-telescope, formerly known as the Karoo Array Telescope (KAT), 
and has also been shortlisted to host SKA [13]. It is foreseen that the meerKAT will be the 
largest radio-telescope in the world, until the SKA is commissioned, thus South Africa will have 
a strategic advance in, amongst others, the study of pulsars and neutron stars. 

KAT-7, a seven-dish engineering test bed, will be commissioned towards the end of 2009 and it 
is foreseen that the full fifty-dish array will be completed by 2012 |13| . 

The development of the meerKAT could not only be of importance to radio-astronomers but 
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also to nuclear physicists studying dense matter. This work can be seen as one example of how 
nuclear physics can be advanced by using data gathered from outer space. 

1.2 Conventions and notations 

In this thesis the following conventions and notations will be used. 

1.2.1 Vectors, fields and operators 

Three dimensional vectors will be expressed in bold font, e.g. k, or using Latin indices, k l . 

Four dimensional vectors will be denoted by the Lorentz (Greek) indices, e.g. fi. In the case of 
the position vector, x^, it is given by 



In most instances the Lorentz index in the four-vector will be suppressed, after it has been 
initially defined with the index. 

The explicit dependence of a field on spacetime coordinates will be suppressed, especially in 
later chapters, i.e. 



If the field or the vector is constant it will be defined as such and denoted by a subscript, usually 
zero e.g. </>o- 

Derivatives with regards to contravariant, x^, or covariant vectors, defined as 



(t> l x{x l/ ) = <?V(z) = </> M • 








d d 
dx° ' dx 1 ' 



d d 

dx 2 ' dx 3 




and 




where the properties of the metric tensor (|2.3|) have been used. 



The time-derivative is sometimes expressed as 



at ~ 
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1.2.2 Natural units 

In this work the notion of "natural units" , where 

h = c = 1, 

will be adopted. This will be used since some equations are very cumbersome if all the factors 
of h and c are included. Where the inclusion of h or c is of specific interest or meaning they will 
still be initially included. 



This convention has the added benefit that it establishes a simple relation between mass and 
length scale, through the relation 

hc = l = 197.33 MeVfm. 



In natural units the relativistic energy-mass relationship, 



E 2 = p 2 c 2 + m 2 c 4 , 



reduces to 



E 2 =p 2 + m 2 , 



and the units of mass and energy are the same. Table [LT1 contains the constants and conversion 
factors (in natural units) used in this work, that are not explicitly defined elsewhere. 

1.3 Matrices 

The N x N identity matrix will be denoted by 1. If the value of N is known or specific, it 
will be denoted by a subscript, i.e. for N = 2, 



1 
1 
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Name 


Symbol 


Value 


Solar mass 


M Q 


1.98892 x 10 3U kg or 






1.1155 x 10 60 MeV 


Gravitational constant 


G 


6.6726 x 1CT 11 m 3 kg -1 s -2 or 






1.325 x 10~ 42 fm/MeV 


Length 


1 fm 


1 x 1(T 15 m 


Energy 


1 MeV 




Energy density 


1 MeV/fm 3 


1.7827 x 10 12 g/cm 3 


Pressure 


1 MeV/fm 3 


1.6022 x 10 33 dyne/cm 2 



Table 1.1: Constants used in this thesis. 



CHAPTER 2 
Relativity 



Due to the central role that gravity plays in binding neutron stars, any description of these 
objects must include the effects of gravity. In this work gravitational effects will be included by 
using a description that is consistent with the special and general theories of relativity. 

2.1 Introduction 

One of Albert Einstein's phenomenal achievements was to recognise that mass and energy are 
related. As such both mass and energy are gravitational sources and not only mass, as is implied 
in Newtonian theory. 

Since neutron stars are large, ultra-dense objects with large gravitational fields and high particle 
energies, relativistic effects need to be considered in the description of these objects. The aim of 
this chapter is to give a brief introduction of the concepts and the necessary tools in the general 
and special theories of relativity that are applied in the description of neutron stars. 

2.2 Special theory of relativity 

The principle of relativity states that the same mechanical laws apply in any frame of 
reference. This principle was already clearly formulated by Galileo Galilei in the 1600's for 
classical mechanics (i.e. mechanics governed by Newton's Laws), but Albert Einstein, in the 
special theory of relativity, extended this principle to include all physical laws. [H] 

The special theory of relativity (commonly referred to as special relativity) is grounded upon 
two postulates [T5] : 

The postulate of relativity: That the same physical laws apply in any inertial (uniformly 
moving) frame of reference. 

The speed of light is universal: The speed of light in a vacuum (denoted by c) is the same 
for all inertial observers, regardless of the motion of the source. 

The line element in special relativity, known as the proper time dr, is given by 



dr 2 = r]^ v dx^dx v = c 2 dt 2 — dx 2 — dy 2 — dz 2 . 

8 
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It is the preserved interval that defines Minkowski space, refers to a spacetime point and is 
expressed in Cartesian coordinates as 

x» = (ct,x,y,z). (2.2) 
From the expression of the line element (|2.ip . the metric tensor of Minkowski space is given by 



1 














-1 














-1 














-1 



(2.3) 



Special relativity does not include any gravitational effects and therefore spacetime is flat and the 
metric tensor (|2.3j) is fixed. Frames of reference in which the expressions of special relativity are 
valid are called Lorentz (inertial) frames. These frames are related by Lorentz transformations. 
Due to the postulate of relativity the physical expression of any observable should not differ 
when observed from different Lorentz frames. Because any two inertial frames of reference are 
related by a Lorentz transformation it means that the equations describing a physical observable 
must be invariant under Lorentz transformations. This property is known as covariance and if 
an equation is expressed in a covariant manner the form of the equation does not change when 
transformed from one Lorentz frame to another. To be able to write equations in a covariant 
manner tensor analyses is used. [14] 

2.3 General theory of relativity 

The special theory of relativity gives a very accurate description of our world within a Lorentz 
frame, but neglects one of the most prominent forces that govern the large scale behaviour in 
our universe, namely gravity The distinction between gravity and other forces is that gravity 
acts on all matter, regardless of their internal structure or composition: Particles with wholly 
different physical properties will follow the same free path in a gravitational field, given only 
that their initial velocities are the same. The general theory of relativity (commonly referred to 
as simply general relativity) expands the locally very accurate special theory of relativity to a 
global scale, by including gravitational effects. [2] 

Special relativity would unfortunately not suffice if gravitational effects are included since no 
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global Lorentz frame can be constructed in a non-uniform gravitational field |14j . This is because 
Euclid's Parallel Postulate holds in Minkowski space, i.e. in Minkowski space parallel worldlines 
stay parallel if extended infinitely. In a non-uniform gravitational field particles do not travel 
along straight lines, but along curved trajectories: if two particles are released from the same 
height above the earth's surface, but some distance apart, they will both fall towards the centre 
of the earth. Their trajectories would be parallel on a local scale but overall they will follow a 
curved path [4]. The spacetime on a global scale and in the presence of gravitational sources is 
one in which Euclid's Parallel Postulate does not hold, i.e. a spacetime that is flat on a local 
scale but overall curved, as described by Riemannian geometry |14| . 

The absence of a global Lorentz frame might spell trouble for special relativity: if no Lorentz 
frame with regards to an event can be found, special relativity will be a very interesting but 
unphysical theory. Luckily within a free-falling frame in a non-uniform gravitational field the 
effect of the gravitational field would not be noted. Thus a free-falling frame of reference would 
constitute a Lorentz frame with regards to the event. Einstein used this fact to relate the special 
and general theories of relativity through the Equivalence Principle: 

The Equivalence Principle: Physical laws that describe how the forces of nature behave 
in a gravitational field have the same (covariant) form as in the special theory of relativity when 
these laws are expressed in a frame of reference, free-falling with the event, in the gravitational 
field. 

The Equivalence Principle is analogous to the theorem in differential geometry that states that 
a flat space tangent to the manifold can be imposed on any point in a differentiable manifold. 
The spacetime of general relativity constitutes such a manifold and therefore a Lorentz frame 
can be imposed at any point. [H] 

To translate from special to general relativity the metric tensor, rf v , has to be replaced by 
the general metric, g^ u , and normal derivatives by covariant derivates (see Sec. 12.41 for details of 
g^ u and covariant differentiation) [4]. 

2.4 Riemannian spaces 

Albert Einstein recognised the similarities between gravitational physics and Riemannian 
spaces and used the mathematics describing these spaces to describe gravitational phyics. 
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In this section the main points of Riemannian spaces will be introduced. 
2.4.1 Geodesies 

A line in curved space is defined as straight (i.e. the shortest distance between two points) if 
the tangent vector to the line is parallel-transported [2]. These lines are known as geodesies. 



The metric tensor of curved space is denoted by g^ v and is defined in terms of the natural 
basis vector (g ) of the space 

ds 

9 



where s is a small displacement along the coordinate curve. The metric tensor is given by the 
dot product of two basis vectors 

In Minkowski space the metric tensor (rjuv) is fixed and given by (|2.3|) . In curved space the 
metric tensor is not necessarily constant, but dependant on the geometry of the space. [IT] 

Geodesies, just as straight lines in a Euclidean space, can also be defined as the extreme value of 
the integral of the proper time interval (dr) between the two spacetime points. The line element, 
the proper time, in general relativity is given by 

dr 2 = g^dx^dx". (2.4) 

As is shown by N.K. Glendenning in Ref. [I], geodesies are thus described by the Geodesic 
equation, 

d 2 x x A dx»dx" 

where is the affine connection. It also shown in Ref. [1] that the affine connection is equal to 
the Christoffel symbol of the second kind. 



2. Relativity 



12 



2.4.2 Christoffel symbols 

The Christoffel symbol (of the second kind), T" , is defined in terms of the relation of the 
derivative of the natural base vector to the natural base vectors as |17j . 

M 

As shown in Ref. [16j . the Christoffel symbols can be expressed in terms of the metric tensor 

and so it can be seen that the Christoffel symbols are related to the curvature of the space 
through the change in the metric. 

The Christoffel symbol of the first kind, T KUfJi , is related to the Christoffel symbol of the second 
kind through [5J: 

The Christoffel symbols of the second kind are symmetric in their covariant (lower) components 
[IT], i-e. 

■pet -pa 

and so the Christoffel symbols of the first kind are symmetric in the last two covariant compo- 
nents. 

2.4.3 Covariant derivatives 

The partial derivative of a vector v, with v = v^g , is: 



dx v 

v %g„ + g„ „ v 



= —g + 
dx v dx u —^ dx v 



p, —p,y 

= <^ + ^r^ from (USD 

and the covariant derivative is defined as [TTj: 

< = < + ?/T^. (2.7) 



2. Relativity 



13 



The covariant derivative transforms covariantly, i.e. the covariant derivative of a tensor is once 
again a tensor [16] . 

2.4.4 Riemann-Christoffel curvature tensor 

The Riemann-Christoffel curvature tensor (or just curvature tensor) describes the curvature 
of space. It is given by 

OP _ pP _ pP I p« pP _ pa pp fry o\ 

-""cr/W cti/jP 1 erf o^t - 1 - trp- 1 - cm/" V^ ,l -V 

The curvature tensor is derived by evaluating the parallel-transport of a tangent vector along a 
small closed loop on a curved surface. The vector field is defined on the surface and a vector 
(V**) parallel-transported along the loop. Once the starting point is reached again, the vector is 
compared to the original vector. The difference between the two vectors is proportional to the 
curvature and the curvature tensor |14| . 

A more intuitive explanation might be given when it is considered that in flat spacetime the 
metric tensor is constant and the order of covariant differentiation does not matter. But in 
curved spacetime the order of covariant differentiation is important as two successive covariant 
differentiations do not commute, as is shown in Ref. (141 : 

where the commutator is given by the curvature tensor. The commutator is analogous to com- 
puting the change in a vector along a closed path: the change in the vector is first computed in 
one direction and then in another and the changes in the reverse order is subtracted. |14| 

The curvature tensor is zero if and only if the space is flat |16| . 

2.4.5 Bianchi identities 

From the properties of the curvature tensor it is shown in Ref. [H] that the following equation 
holds in any frame: 

Ra/3fii/;X + Ra/3\n;u + Ra(3v\;[i = . (2.9) 

This is known as the Bianchi identities. These identities can be used to determine the metric 
tensor if the curvature is known [16J. 
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2.4.6 Ricci tensor and scalar curvature 

The Ricci tensor, R^ u , is a contraction of the curvature tensor: 

Ra/3 = Rav(3 = Rpa- (2.10) 

The Ricci tensor can also be formed by contracting indices, other than the first and the third, 
but the symmetry of the curvature tensor implies that these contractions will either vanish or at 
most add a minus sign to the Ricci tensor. The Ricci tensor is symmetric. 



The scalar curvature, R, is defined as the contraction of the metric tensor and the Ricci tensor 

R = g^R^v = g^R^w 

2.4.7 Einstein tensor 

If the Bianchi identity is contracted twice |14| : 

9^ [RaPiiv;\ + Ra/3Xfj,;v + Ra(3uX;^} = R/3u;X ~ R(3X;i/ + Rp v X;n = ® 

and 

9 R/3u;X ~ R/3X;u + Rp v X;n = R \ x ~ R ^ ~ R X^ = °' 

then the last equation can be rewritten as: 



By defining as 



(2^-^i?). M = 0. (2.11) 



GT = jff» - \^f v R (2.12) 



it is clear from Eq. (|2.11[) that it has vanishing covariant divergence, i.e. 

= • (2.13) 

G^ u is known as the Einstein tensor and it plays a fundamental role in general relativity (see 
Sec. 12. 6p . From the symmetry of the Ricci and metric tensors it is clear that the Einstein tensor 
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is also symmetric. 

2.5 Energy-momentum tensor 

The energy-momentum tensor (T^ u ) describes the internal properties of an energy-mass 
distribution (i.e. matter) and is the source of the curvature of spacetime. To be able to give an 
accurate, covariant description of the internal properties of matter within a relativistic frame- 
work, not only the energy distribution, but also the distribution of momentum within the matter 
needs to be considered. 

The energy-momentum tensor is defined in terms of the flux of the 4-momentum across a con- 
stant surface. The element, T^ u , can be described as the /^-component of the momentum flux 
across the x^-surface. Thus the following components can be explicitly named: 

T 00 = energy density, 

T 0i = energy flux across the x l surface, 

T*° = i momentum density, 

T tJ = flux of the i momentum across the j surface. 

If we evaluate the energy-momentum tensor in a momentarily comoving frame of reference all 
elements of the matter would be momentarily static (no spatial momentum), but that would not 
imply that there is no transfer of energy. Energy might be transferred by heat conduction and 
therefore the T°* terms of the energy-momentum tensor would be non-zero. The T*° components 
would also be non-zero for the energy that is being transferred will carry momentum. The T° l 
components would be equal to the T i0 components since the energy flux is the energy density 
times the speeds at which it flows. In the relativistic frame, mass is equal to energy so the 
energy flux is equal to the mass density times the speeds at which it flows, which is the density 
of momentum. 

The spatial components of the energy-momentum tensor (T y ) are symmetric. If they were not 
it would be imply that the elements are whirling around inside the fluid, in the momentarily 
comoving frame of reference. Thus T^ v is in general symmetric and symmetry of the tensor in 
one frame of reference would imply that it is symmetric in all frames of reference. 
The spatial components of the energy- momentum tensor represent the forces between adjacent 
elements in the matter: the off-diagonal spatial elements represent the viscosity (or any other 
forces that are parallel to the interface between elements) and the diagonal elements, the normal 
forces (i.e. the pressure). [14] 
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In general the energy-momentum tensor for static, spherically symmetric perfect fluid (no vis- 
cosity or heat conduction) , moving with a velocity v is [U [21] : 



T» v = -Pjf v + {P + e)u^u u , (2.14) 



where 

• e is the energy density, 

• P is the pressure, and, 

• vF is the four-velocity, 



dx* 



Vl-v 2 (l,v\v 3 ,v 3 ) , (2.15) 



and therefore u^u^ = 1 [3]. 



The conservation of mass and energy is expressed in terms of the energy-momentum tensor by 
the fact that the energy-momentum tensor has vanishing divergence, i.e. 

T^ :K = 0. (2.16) 

In general relativity Eq. (|2.16p generalises to [16] 

T mK = 0. (2.17) 

2.6 Einstein's field equations 

Physical theories are mostly categorised by defining equations, which in the case of general 
relativity are the Einstein field equations. The Einstein field equations are given by solutions 
of the Einstein tensor. For spacetime inside a distribution of mass and energy (i.e. a star), 
Einstein's field equations are given by 

CT = kT^, (2.18) 

where T^ v is the energy-momentum tensor and k is a constant that is determined by comparing 
general relativity to Newtonian mechanics in the Newtonian limit [3j. The energy- momentum 
tensor is a symmetric, divergenceless tensor that is constructed from the mass-energy properties 
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of the medium. 

The Einstein field equations (I2.18P do not only govern spacetime within a mass-energy dis- 
tribution, but also the arrangement of mass-energy distribution as well as its dynamics. Thus 
spacetime is acted upon by a mass-energy distribution, as these are sources of gravity. On the 
other hand, spacetime influences the mass-energy distribution through the curvature of space- 
time. Thus spacetime and matter co-determine the universe within which we live. 



CHAPTER 3 
Tolman-Oppenheimer-Volkoff equations 



By using the concepts defined in Chapter [2] a relativistic description for a neutron star will be 
derived in this chapter. 

3.1 Introduction 

As a first approximation a star (such as a neutron star) is assumed to be a static, spherical 
symmetric fluid (gas) in hydrostatic equilibrium [3]: the star is bound by an external pressure 
that compresses the star (gravity) while the star is stabilised against gravitational collapse by 
the internal pressure in the star due to some repulsive force. 

If a neutron star is assumed to be a static spherical symmetric fluid in hydrostatic equilibrium 
it would have to be described by a relativistic equation for hydrostatic equilibrium. 

This problem was first studied by R. C. Tolman [18] and by J. R. Oppenheimer and G. M. 
Volkoff [19J in 1939. They derived an equation to study hydrostatic equilibrium in a relativistic 
environment by describing a neutron star assuming that it consists of a neutron gas at high 
density [7J: 

dPjr) Ge(r)M(r) 
dr c 2 r 2 

where 

dM(r) _ 4ire(r)r 2 
dr c 2 

and 

• M(r) is the enclosed mass of the star, 

• e(r) is the energy density, and 

• P(r) is the internal pressure of the star. 
Eq. (|3.ip is known as the Tolman-Oppenheimer-Volkoff (TOV) equation [20J. 



1 + 



P(r) 
e(r) 



1 + 



A-Kr 6 P{r) 
M(r)c 2 



2GM(r) 
c 2 r 



(3.1) 



The TOV equation expresses how the pressure decreases from the centre of the star to the 
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edge (where P = 0) in terms of the energy density and the pressure. These two quantities are 
related through the equation of state of the matter in the interior of the star. The larger the star, 
the higher the pressure in the centre of the star will be, since it is gravity that compresses the 
star. Since the pressure appears on the right-hand side of Eq. (13. ip . an increase in the central 
pressure will increase the pressure gradient [the left-hand side of Eq. (|3.1j) ]. Thus the more 
massive the star, the smaller the radius at which P = 0, since the pressure gradient is steeper. 
Thus there exists a critical value for the mass of the star, above which the star will collapse 
under its own gravity [4] . This critical value is known as the maximum mass of the star and will 
differ for different equations of state of the matter in the neutron star interior. 

3.2 Derivation of the TOV equation 

The following derivation is taken almost entirely from Ref. [4]. (The derivation in Ref. [4] is 
more elaborate and hints to perform some of the calculations are also given.) 

To derive the TOV equations from the Einstein field equation (I2.18h it suffices to consider a 
static, isotropic star. The line element for static, isotropic spacetime can be given by: 

dr 2 = e 2 ^dt 2 - e 2X ^dr 2 - r 2 d9 2 - r 2 sin 2 Qd<\? (3.3) 

where v{r) and A(r) are functions that need to be determined and x^ is the spacetime point 
described in natural units by 

x" = [t,r,0,<l>]. 



From the expression for the line element (|3.3p the components of the metric tensor can be read 
off as 

500 = e 2 ^ (3.4a) 

5ii = -e 2X(r) (3.4b) 

522 = ~r 2 (3.4c) 

533 = -r 2 sm 2 9 (3.4d) 

<T = <V = (for/i^i/). (3.4e) 
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Since g^ a g au = S", g^y is given by g^ = (g^) 1 (no summation over indices). 

Using the symmetry properties of the Christoffel symbols, the components of the Ricci tensor in 
static isotropic spacetime are: 

9;/ 

R 00 = (V + AV-!^-- )e 2 ("- A \ (3.5a) 

r 

9 V 

R u = I/ "_AV + Z / /2 -— , (3.5b) 

r 

i? 22 = (1 + ru' - rA')e 2(A) + 1, and, (3.5c) 
i? 33 = R 22 sm 2 e, (3.5d) 

where the primed index refers to differentiation with respect to r. 

For the construction of the TOV equation it is easier to work with mixed tensors, therefore 
the Einstein field equations (|2.18p are written as 

= g^ a R au — -^g^ag au R 
= Rf. Sf.R 

= kT;, (3.6) 

and thus, using the properties of the Ricci tensor (|3.5|) . the components of the Einstein tensor 
become: 

r 2 G° Q = e" 2A (l - 2rX') - 1 = -— [r(l - e~ 2X )] (3.7a) 

dr L 

r 2 G\ = e- 2X {l + 2rv')-\ (3.7b) 

/ 7 J \l 

i2 -2A I ,11 , ,/2 \/../ , v A 



G\ = e- 2X \u" + v' 2 - AV + —— - j (3.7c) 
Gl = G\. (3.7d) 

From (12. 14ft and the equivalence principle, the energy-momentum tensor is 



= -P(r)gV + (P(r) + e(r))uV. 
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Since the the star is static (u 1 = 0) and therefore from (|2.15j) . g^u^u" = 1, u° is given by 

1 



u° 



V9oo 

From the equivalence principle, a comoving Lorentz frame can be defined 0, and therefore the 
energy-momentum tensor can be expressed as 

T ° = e(r) and T? = -P(r). (3.8) 



From (jTTajl . (EHJ) and (JM 



dr - 
K r 2 T ° 



'(1 



3 -2A(r). 



Kr 2 e(r) 



(3.9) 



and (|3.9|) can be solved to obtain 



e 



-2A(r) = 1-1 e {r)r 2 dr 



r 



o 



where M(r) is defined as 



1 - £-M(r), (3.10) 



M(r) = 4vr / e{i J )i J2 di J . (3.11) 



o 



Eq. (|3.1ip is an expression for the gravitational mass i.e. the mass-energy that generates a 
gravitational field included, up to radius r, in the star. The relativistic expression (in natural 
units) that is used to describe the relation between the energy density (e) and the mass density 

(p) is: 

e(r) = p{r). 



lr This is a valid assumption, since the change in the metric from the centre of the star to its boundary is not 
significant over the spacing between a few nucleons [4] . [The change in the metric can be estimated by considering 
the espression for gu (|3.14l) .] 
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Using (|3.10p the following expression for (|3.4bj) can be obtained: 

9u = -eW = -(l-"-my\ (3.2) 
For ()3.12p to agree with Newtonian mechanics in the Newtonian limit, k should be defined as 

K = -8ttG, (3.13) 
with G the Gravitational constant [3]. Thus (|3. 12|) becomes 



5ii = -e 2A{r) = - ( 1 - ) ■ (3.14) 



2GM (r) 



-i 



Using (I3.13P and (13.8P the expression for the Einstein field equations (13. 7p . becomes: 

/ i 2A'\ 1 

G°o = e- 2A - —J - - 2 = -8nGe(r) (3.15a) 

G\ = e~* ( i + *£\ - 1 = -8vrGP(r) (3.15b) 



G\ = e~ 2X ( v" + v' 2 - AV + ) = -8vrGP(r) (3.15c) 



C3 = G2 (3.15d) 

By manipulating (I3.15ap . (I3.15bp and (I3.14|) and substituting it into (|3.15cp . the following ex- 
pression for P'{r) is obtained: 

G[e(r) + P(r)] [M(r) + 4^r 3 P(r)] 
^ (rj "dr" r[l-2GM(r)] ' [6Ab) 

Eq. (|3.16p is the TOV equation in natural units. 

3.3 Solving the TOV equation 

To obtain sensible information from the TOV equation the coupled differential equations, Eqs 
(|3.ip and (|3.2p . must be solved. To be able to do this the initial conditions (at r = 0) for P(r) 
and M(r) must be known. 



Since the expression for M(r) ([3. lip refers to the enclosed mass in the star and r = refers 
to the centre of the star, the initial values for both M(r) and r are taken to be zero (or very 
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small). The choice for the initial P depends on the inferred central baryon density of the neutron 
star. 

Using the initial values for P, r and M the pressure (P(r')) at an incremental increase in r, 
namely r' = r + dr can be calculated using Eq. (13. lj) . Using the equation of state, the cor- 
responding energy density, e(r'), can be calculated. Using ()3.2|) the value for M(r') can be 
obtained. This cycle is repeated until r' = R at which the P(r') is zero, is reached. This denotes 
the boundary of the star and R is the radius of the star. M{R) is the mass enclosed by the 
radius of the star, i.e. the mass of the star. 

Each different assumption of the central density will correspond to a unique mass and radius 
relationship of the neutron star. Therefore by varying the initial central density, a whole sequence 
of possible masses and radii of neutron stars, corresponding to a specific equation of state, can 
be generated. 

As was expounded upon at the beginning of this chapter, each neutron star sequence has a 
maximum mass and as such each equation of state of the neutron star interior implies a different 
maximum mass. The equations of state with higher maximum masses are referred as stiff, i.e. 
the pressure increases rapidly with an increase in density. The equations of state with lower 
maximum masses are referred to as softer equations of state [1]. 



CHAPTER 4 
Modelling dense nuclear matter 



The challenge in the description of matter at high densities, such as that in neutron stars, is 
to develop a model that not only describes matter at high densities, but also the properties 
of matter observed at normal densities. Additionally, relativistic effects become much more 
pronounced in dense systems as the particles attain energies that are comparable to their rest 
mass. Therefore it is neccesary that a description of matter at high densities must incorpo- 
rate the general properties of quantum mechanics, Lorentz covariance, electromagnetic gauge 
invariance and microscopic causality within a many-body system [21]. The only framework that 
can describe a relativistic, quantum-mechanical, many-body system in such a way is relativistic 
quantum field theory based on a local, Lorentz-invariant Lagrangian density [21.) 22] . 

Quantum chromodynamics (QCD) describes the interaction between quarks via the exchange 
of gluons. This fundamental theory is an obvious candidate to describe dense matter systems, 
but unfortunately this theory has computational difficulties at such scales [21] . It is also cumber- 
some to describe nuclear phenomena in terms of quarks and gluons since quark degrees of freedom 
are not observed in nuclear experiments but, hadronic degrees of freedom. Hadrons are particles 
that consist of a number of quarks and are further subdivided into baryons and mesons: baryons 
are particles that contain three quarks (such as protons and neutrons), while mesons contain a 
quark - anti-quark pair. A description of the interaction between two nucleons (particles in the 
nucleus), based on the exchange of mesons was first introduced by Hideki Yukawa in 1935 [23] 
and gained much acceptance after the actual discovery of mesons. A complete discussion of the 
mesonic theory of the nucleus is given in Ref. [23] . 

A quantized field theoretical description of nuclei and nuclear matter, called quantum hadrody- 
namics (QHD), which is based on hadronic degrees of freedom was introduced by John Walecka 
in 1974 [23]. It should be noted that QHD is not a fundamental theory, but an effective one since 
hadrons are composite particles. 

In QHD, as in QCD, one soon runs into computational difficulties and to obtain solutions from 
the description certain approximations have to be made. The coupling constants in QHD are 
large and thus a perturbation expansion in terms of the coupling constants will not suffice. In- 
stead the system can be approximated in the relativistic mean-field. This approximation will 
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be discussed in Sec. 18,31 The saturation properties of nuclei and nuclear matter, that any de- 
scription of dense matter should reproduce, are discussed in Section 14.11 In Chapter [5] a brief 
overview of relativistic quantum mechanics is given. This chapter is followed by a brief discussion 
of the necessary tools and concepts in classical field theory that is used in quantum field theory 
in Chapter [H Quantum field theory is discussed in Chapter [7j after which QHD is formulated 
in Chapters M and [H 

4.1 Properties of symmetric nuclear matter constraining nuclear models 

Symmetric nuclear matter, or just nuclear matter, is an idealised system that stems from one 
of the original models of the nucleus, the liquid-drop model. The properties of nuclear matter 
are inferred from the experimentally observed properties of finite nuclei [3J, but, since it cannot 
be directly observed, there seem to be some disagreement as to what the exact values of certain 
properties of nuclear matter are. In this work the values of the properties of nuclear matter given 
in Ref. [4] will be taken as the observed values. 

4.1.1 Saturation density 

Nuclear matter is a saturated system, due the characteristics of the strong interaction. The 
short-ranged, strong nuclear interaction is the dominant interaction between nucleons. It is 
essentially attractive, which is necessary to form stable nuclei, but repulsive at short distance 
(< 0.4 fm). Since the strong force acts only over a short distance, the interaction is limited to 
nearest neighbours in a dense system. Therefore at a certain density the central density of the 
system will not increase any further, even as more nucleons are added to the system. The density 
at which this occurs is referred to as the saturation density. At saturation density the pressure 
of the system is zero and the system will remain in this state if left undisturbed. [4J 

The density of saturated nuclear matter given as 0.153 fm~ 3 in Ref. @j and 0.16 fm" 3 in Ref. [6]. 

4.1.2 Binding energy 

In general terms the binding energy of a system is the energy expended or required to form a 
system. In the case of a stable system the binding energy is negative, hence energy was donated 
to the surrounding environment in forming the system, leaving the system at a lower energy 
state than the sum of all the individual parts of the system. 

At saturation density the binding energy of the system will be at a minimum since at this 
density the system is in its most stable (lowest energy) state, compared to other densities. The 
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binding energy of nuclear matter is given as -16.3 MeV/nucleon in Ref. [4] and -16.0 MeV/nucleon 
in Ref. [6] 

4.1.3 Symmetry energy 

Stable nuclei with low proton {Z) number prefer a nearly equivalent neutron (N) number. 
As Z increases the (repulsive) Coulomb interaction between the protons also increases. Stable 
nuclei then diverge from being N = Z nuclei to nuclei with a N larger than Z, as can been seen 
on a Segre diagram. The influence of this preference is accounted for by the symmetry energy. 
The symmetry energy coefficient, 04, stems from liquid-drop model of the nucleus, and refers to 
the contribution made by the isospin assymetry to the energy of the nucleus [25J. 

The symmetry energy coefficient is given by: 



The value of 04 is estimated to be between 31 and 33 MeV according to Ref. [5], while Refs [3] 
and [6] state the value of 04 to be 32.5 MeV (without any uncertainty). 

4.1.4 Compression modulus 

The compression modulus defines the curvature of the equation of state at saturation [4J and 
is related to the high density behaviour of the equation of state. If the energy density rapidly 
increases with an increase in pressure the equation of state is referred to as stiff. With a soft 
equation of state the energy density increases more gradually with an increase in the pressure 



The value of K has been estimated to 234 MeV (with some uncertainty) 0j. Ref. [6] state the 




(4.1) 





(4.2) 



value of K to be 265 MeV. 



CHAPTER 5 
Relativistic quantum mechanics 



This chapter is a very concise introduction to relativistic quantum mechanics, introducing only 
concepts necessary within the scope of this work. A much more thorough discussion can be found 
in Refs [26] to [30]. 

5.1 Introduction 

Relativistic quantum mechanics is the merger of quantum mechanics and special relativity. 
One possible candidate for a relativistic quantum mechanical description could be the Schrodinger 
equation. The Schrodinger equation can be motivated by applying the canonical quantization to 
the non-relativistic energy-momentum relationship |26j . 

2m 

where the physical quantities are replaced by operators, i.e. 

d 

E -> ih— (5.1a) 
dt K ' 

p -> -ihV, (5.1b) 

and by letting the operators operate on the wave function ijj(t,x), the Schrodinger equation is 
obtained: 

.a . -h 2 v 2 ,, . 
ih—^{t,x) = ——^(t,x). 



Relativity requires that the equations of motion in one inertial frame are valid in all inertial 
frames, meaning that all equations describing the motion of a particle must transform covariantly 
under Lorentz transformations. The Schrodinger equation does not transform covariantly and 
therefore another approach is needed to describe relativistic particles in a quantized theory |27j . 
An obvious candidate would be to apply the canonical quantization to the relativistic mass-energy 
relationship. 
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5.2 Klein- Gordon equation 

Imposing the canonical quantization procedure to the relativistic energy-momentum relation- 
ship, 

E 2 = p 2 +m 2 , 



results in the Klein-Gordon equation 

-|^) = (-V 2 + m 2 )<K^), 
where 4>(x*) is some wave function that depends on the four-vector x* = (t,x). 

The Klein-Gordon equation can be expressed in a covariant manner as: 

(dud* + m 2 ) <f){x u ) = 0. (5.2) 



Since 

^ ~ an 

transforms like a scalar under Lorentz transformations [3], the wave function, <j), has to describe 



8 2 

duff* = ^- - V 2 



a scalar particle for Eq. (|5.2p to be a covariant equation. Therefore plane-wave solutions of the 
form 

oc e - iEt+ip x = e~ ip - x , (5.3) 

with p^ 1 = (E,p) the momentum four-vector and p ■ x = p^x* = Et — p ■ x, can be solutions of 
the free particle Klein-Gordon equation [261. 

Substituting this wave function into the Klein-Gordon equation (|5.2p . the relativistic energy- 
momentum relation, 

E 2 =p 2 + m 2 , 
is obtained. As such, the possible energies are: 



E = ±^p 2 + m 2 . 



Thus the Klein-Gordon equation allows for positive as well as negative energy solutions. 
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The Klein-Gordon equation has the conserved current j^ 1 (i.e. d^j^ = 0), with = (p,j), 
where 

p = i{4>*d 4> - 0<9 o </>*) 

j = _i(0* V( />- W)- 

If the time-like component of the conserved current (p) is interpreted as the probability density, 
the negative energy solution would imply that the probability density will not be positive definite 



Historically the possibility of negative energy solutions and the probability not being positive 
definite led to the Klein-Gordon equation being abandoned as a physical description [26]. How- 
ever, the Klein-Gordon equation was redeemed when a physical interpretation was given to the 
negative energy solutions (see Sec. [51 



5.3 Dirac equation 

P. A. M. Dirac wanted to formulate an equation based on the relativistic energy-momentum 
relationship that is linear in spatial and time derivatives in an effort to overcome the difficulties 
of the Klein-Gordon equation, especially the positive definiteness of the probability density 



One way to obtain the Dirac equation, using the relativistic energy-momentum relationship, 
is by making the following substitution [29 : 



E = \/ p 2 + m 2 — > a ■ p + /3m, 

where a and (3 need to be determined, but are assumed to be constants, independent of spacetime 
and to commute with position and momentum operators [30j. 

Making the substitutions of the canonical quantization (|5.ip . the above expression reduces to 
the Dirac-equation 

= (-ia • V + /9m)^(a?) . (5.4) 

As is discussed in Ref. [27) the properties of the Dirac equation (|5.4p can be derived by considering 
certain requirements: Since Eq. (|5.4p must be invariant under spatial rotations, a and (3 cannot 
simply be numbers. If it is considered that the probability density, ip*ip, must be the time 
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component of a conserved four- vector if integrated over all space at constant time, Dirac proposed 
that the wave function must be analogous to a ^-component spin wave function, 



(5.5) 



The column vectors satisfying the Dirac equation are called Dirac spinors. Therefore a and (3 
must be N x N matrices. Since the Dirac equation ()5.4|) must satisfy the relativistic energy- 
momentum relationship, the following properties can be derived: 



{a;,/3} = i = 1,2,3 
{cti,aj} = 28ijl i,j = 1,2,3. 



(5.6a) 
(5.6b) 



Since the Hamilton operator, in this case 



a ■ p + (3m, 



must be hermitian it implies that a and f3 must also be hermitian [30J. Thus a and j3 both must 
be traceless, even dimensional matrices with eigenvalues of dbl. 

For N = 2 the conditions of (j5.6[) are satisfied by the Pauli-matrices and the identity matrix, but 
not all these matrices are traceless. Thus the lowest dimensionality that satisfies all conditions 
for a and (3 is N = 4 \27\ [28] . One representation of these matrices is known is the Pauli-Dirac 
representation, given by : 



<Ti 

a*, o 



12 







-1 2 



10 
10 
0-10 
-1 



(5.7a) 



(5.7b) 
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where 



(7 2 



0"3 



1 

1 

-i 

1 

1 

-1 



, and, 



(5.8a) 
(5.8b) 
(5.8c) 



are the Pauli-matrices. 



In the non-relativistic limit the Dirac equation reduces to the Pauli equation describing spin-^ 
|27j . Thus the free particle solutions of the Dirac equation can be found in analogy to that of a 
non-relativistic spin-^ particle, where the wave functions consist of a two component spinor and 
a plane- wave function [26]. A possible solution to the Dirac equation (j5.4|) is [26J 



t/)(x) = ujc tp '- 



(5.9) 



where and are spacetime and momentum four-vectors and w is a four component column 
vector written as 



X 



with (j) and x two component spinors. 
Solving for ip, Eq. (|5.4p reduces to: 



E \ 


( 4>\ 











m\-2 cr • p 
<j ■ p —ml2 




thus obtaining coupled equations for (j) and x : 



(E — m) (ft = cr ■ px 
{E + m) x = <r ■ p<f> . 



(5.10a) 
(5.10b) 
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From Eq. (|5,10b|h oj is given by: 



to 



E+m 



When the solution for \ [from Eq. (|5.10b|) ] is substituted into Eq. (|5.1Uap . it can be shown that 
[2BJ: 

(E - m)(E + m)4> = p 2 (f>, 

which implies that 

E = ±\/p 2 + m2 i 

thus allowing positive and negative energy solutions to the Dirac equation [26j . 

To be able to write the Dirac equation in a covariant form, the 7-matrices are defined as: 

7 = p (5.11a) 
7 < = ptf = 7°a*. (5.11b) 

The properties of the 7-matrices in the Pauli-Dirac representation are derived from the properties 
of a and /3 (|5.6p and that of the Pauli-spin matrices f|5.8j) : 



{7^,7"} = 2^1 4 , 

(f) 2 = -1 4 for i = 1,2,3. 

(7°) 2 = 14, 



and the adjoints of the 7-matrices are 



y4 _ 7 ° 7 ^ 7 °. 



The Dirac equation (|5.4p can be written in a covariant form, by multiplying it from the left by 



7 , as: 



(ry"<9 M -m)il>{x) = 0. (5.13) 
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Eq. (|5.13p will only be covariant under a Lorentz transformation, 



x* 1 A (x u y = a;x , 

if the wave function, if)(x) transforms as 

4>(x) ^ ^'{x') = S{A)4>(x) , 

where S(A) is some matrix |28j . 

As shown in Refs [27] and [28], S(A) must satisfy: 

5- 1 (A) 7 ^(A)= 7 -. 

Defining the Dirac adjoint ip as 

^ = ^(x)-f°, (5.14) 
it is also shown in Ref. [27] that tp transforms under a Lorentz transformation as 

Tp'(x') =ip(x)S~ 1 (A). 

The Dirac equation has a conserved current, j^, 

of which the time-like component, p, is given by 

p = Tp^ip = ifj^tfj 

so that 

4 
i=l 

Thus the probability density of the Dirac equation is positive definite, but negative energy 
solutions are still allowed 
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5.4 Interpretation of negative energy solutions of the Klein-Gordon and Dirac 



There are two interpretations of the negative energy solutions. First P. A. M. Dirac put 
forward an interpretation of the negative energy solution of the Dirac equation for fermions and 
then R. P. Feynman added another interpretation that could be applied to fermions as well as 
bosons (scalar particles, described by the Klein-Gordon equation) [29] . 

5.4.1 Dirac's interpretation 

The Dirac description states that there are positive and negative energy states available to 
a free spin-| particle, such as an electron. These energy levels are symmetric around the zero 
energy state. To prevent positive energy particles from spontaneously decaying to negative energy 
states Dirac postulated that in the vacuum state all the negative energy states are filled (the 
Dirac sea [261). Since all negative energy states are filled, Pauli's Exclusion Principle prevents 
the decay of a positive energy electron to a negative energy state. [26] 

An electron occupying a negative energy state can be excited to a positive energy state, leaving 
a hole in the Dirac vacuum. Due to the absence of a negatively charged electron the hole will 
behave as a positively charged particle with a positive energy with regards to the Fermi sea, 
predicting the existence of a positron [26J. 




It should be noted that due to this interpretation the Dirac equation no longer describes a 
single particle state, but a many-particle state, thus necessitating a many-particle description 
such as quantum field theory to describe relativistic particles [27] . 

5.4.2 Feynman's interpretation 

R.P. Feynman interpreted the negative energy state as positive energy particle propagating 
backwards in time or as an anti-particle propagating forward in time [29j. 

The Feynman's interpretation vindicated the Klein-Gordon equation. Since the Klein-Gordon 
equation describes bosons, which do not obey the Pauli Exclusion Principle, the Dirac inter- 
pretation could not be applied to explain the negative energy solutions of the Klein-Gordon 
equation. 



equations 




CHAPTER 6 
Lagrangian mechanics and field theory 



6.1 Lagrangian mechanics 

H. Goldstein et ah explain Lagrangian mechanics in Ref. [31] as a description of a physical sys- 
tem in terms of the degrees of freedom of the system. Lagrangian mechanics does not contain any 
new physics, but is simply an alternative expression of the physical laws governing the equations 
of motion of objects. Since the description is in terms of the degrees of freedom of the system, the 
formulation is not restricted to just describing classical particles, but can also describe discrete 
or continuous systems in the classical or relativistic regimes. Lagrangian mechanics is a desirable 
description since all the information pertaining to the system is contained in one function: the 
Lagrangian. 

The Lagrangian (L) is defined by 



• T is the kinetic energy of the system and 

• V the potential energy of the system, 

and can thus be seen as an expression of the energy of a system. Each degree of the N degrees 
of freedom of the system is described by a generalised coordinate q a , with a = 1, 2, . . . , N. The 
Lagrangian is given in terms of the generalised independent coordinates and the time-derivates 
of these, q a and time, t. 

The generalised momentum (conjugate momentum) of each generalised coordinate is p a with 



The evolution of the system is governed by Hamilton's principle. It states that the system will 
evolve over a time interval [ii,£2] hi such a way that the action (/) is stationary. The action is 
the line integral of the Lagrangian between time ti and time t2- If the action is stationary the 



L = T -V 



where 



dL 
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variation of the action will be equal to zero, i.e. 



= 51 = 5 



L(q,q,t)dt 



The equation of motion of a specific degree of freedom described by the coordinate, q a , is given 
by the Euler-Lagrange equation [5T] : 



dt \dq a J dq a 

The Euler-Lagrange equation follows from the application of Hamilton's principle with regards 
to a specific coordinate |31| . 

6.2 Hamiltonian formulation 

The Hamiltonian formulation is an alternative way to express Lagrangian mechanics. The 
Hamiltonian formulation has found various extension in physics and is therefore also briefly dis- 
cussed. 

The essential difference between the Lagrangian and the Hamiltonian formulation is that in the 
former the Lagrangian is expressed in terms of the generalised coordinates, the time-derivatives 
of the generalised coordinates and time itself, (q,q,t). In the Hamiltonian formulation a new 
quantity, the Hamiltonian, H, is defined and expressed in terms of the generalised coordinates, 
the conjugate momenta of the generalised coordinates and time, (q,p, t). Thus on a mathematical 
level, the switch from the Lagrangian to the Hamiltonian formulation can be seen as a Legendre 
transformation |31j . 

The Hamiltonian is defined by (with the sum over a implied) [31|. 132] : 




H(q,p,t) = q a p a - L(q,q, 



t). 




Anagolous to the interpretation of L, H can be seen as an expression of the total energy of the 
system. If L is L = T — V and if the momentum times the speed of each degree of freedom 
summed over the whole system is twice the kinetic energy (T) of the system, H is given as 
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6.3 Field theory 

In Ref. (31 J a field is denned as set of one or more functions of space and time. These func- 
tions can be used to describe the displacement and velocity of points in a (continuous) system 
|32j . The formalism of previous section can be applied to describe continuous systems using 
fields, by replacing the discrete coordinates, q a , by fields, (f) a . Accordingly the fields 4> a with 
a = 1, 2, . . . , N, describe the degrees of freedom of the system. 

In general, all the fields described in this work will be in 4-dimensional Minkowski space [see Eq. 
(|2.ip for the metric of this space] and thus the fields will depend on spatial and time coordinates, 
i.e. 



In field theory the Lagrangian would be defined in terms of the Lagrangian density (£): 



The Euler-Lagrange equations can be generalised (see for instance Refs |31| and |32| ) to directly 
apply to the Lagrangian density, i.e. 



4> a (t,x,y,z) = (p a (t,x) = 4> a (x^). 



L = 




(6.3) 




(6.4) 



For the field 4> a the conjugate momentum of (|6.ip is defined as TT a |31j 



.a 



dC 



(6.5) 



7T 



•a 
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The Hamiltonian is defined as: 



% dx 
dC 



C ) dx 



= f (fi> a -C)dx, (6.6) 
and thus the Hamiltonian density H is given by 

U = ir a ^ a -C. (6.7) 

6.4 Noether's theorem 

Noether's theorem describes the relation between the invariance of the Lagrangian density 
under certain transformations and conserved quantities (currents). 

The transformation could be a coordinate transformation, 

x » _». x >» = x » + 6x» (6.8) 

or a field transformation, 

<f> a {x) -> <P' a (x') = <t> a {x) + 6<t> a (x). (6.9) 

The theorem states that 

• if the space is flat, and, 

• if the coordinate/field transformation is continuous and the transformation can be written 
in terms of the coordinates as 

SxP = e r Xf, (6.10) 

and in terms of the fields as 

5(j) a = e r ^ ra , (6.11) 

where the functions Xr and ^> rp depend on the coordinate/field variables and e is an 
infinitesimal parameter (r runs from 1,2,...,R), and, 
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• if the Lagrangian density is form- and scale invariant under these transformations, 
then the following will hold |31j : 



I \ a &c" / a si 5 



so that there are R conserved currents of the form: 



dC \ <9£ 

~ 9 CT a (x) - £^ - — — — * ar . 



Form invariance of the Lagrangian density under a transformation implies that 

£ {<P' a (x'),d u cP' a (x'),x') = C (cf/ a (x , ),d v cf> , a (x , ),x') , 
so that the effect of the transformation on the Lagrangian density would be 
C(<j) a (x),d v <j) a {x),x) -> C (4>' a (x'),d u (j)' a (x'),x') . 



Scale invariance implies that the magnitude of the action integral is also assumed to be invariant 
under the transformation, i.e. 



I' = f (dx' 4 )C^' a (x'^,d u( f>' a (x^),x^) 



A proof of Noether's theorem is given in amongst other Refs |31] and |35j . 

6.4.1 Energy-momentum tensor as conserved quantity 

A typical symmetry of the Lagrangian density is translational invariance. This means that 
the Lagrangian density describing the system is independent of the spatial and temporal origin 
of the system. With Noether's theorem one can show that translational invariance imply that 
the energy-momentum tensor is a conserved quantity (also see Sec. 12.5ft [H 1314 [35] . 
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Consider the following spacetime translation 



where hp is an arbitrary constant and e r as in (|6.10|) . Since the Lagrangian density (and therefore 
the fields) are translational invariant, the coordinate change would have no effect on the fields 
and 5(p u (x^) of Eq. (|6.9p would be zero and therefore also ^ ru of Eq. (|6.1ip . 



From Noether's theorem f|6. 12j) it therefore follows that 

d {( dC 



dx u \ \ d{d u 



(PM -C8 v a \h 



Since h a is arbitrary, it follows that 

and so the energy-momentum tensor, T^ u , is given by 

dC 

T ^ = a, x a s ^K-CrT • (6.13) 

In Chapters [8] and [9] the fields, 4> a , that will be considered will actually be field operators. 
Therefore the (ground state) expectation value have to be taken |21j : 

/ dC 
= -Pjf v + (P + e)«V, 

if a static, spherical symmetric fluid, moving with velocity v is consider [from Eq. (|2.14|) ]. If 
v = 0, then 

e = (T 00 ) , (6.14a) 

P = ~ (T u ) . (6.14b) 



CHAPTER 7 
Quantum field theory 



Quantum field theory is the application of quantum mechanics to fields. One of the main mo- 
tivations for this description is the single particle description of the very small scales (quantum 
mechanical) at high (relativistic) energies gives rise to negative energies, which in the interpre- 
tation of Sec. 15.41 gives rise to many-body states. 

Thus to describe this problem in terms of fields has some distinct advantages. The first be- 
ing that field theory is a many-body description in which various degrees of freedom are handled 
with ease. The Lagrangian density is a Lorentz scalar and therefore invariant under any Lorentz 
boosts or rotations. Since the extreme value of the action is therefore also invariant, the Euler- 
Lagrange equation gives rise to Lorentz invariant equations of motion of the fields contained in 
the Lagrangian. A many-particle theory is also necessary in order to preserve causality as is 
shown in Ref. |33j . 

The quantization procedure is very similar to that of normal quantum mechanics, except that in 
this case it is not the particles that are quantized, but fields describing the particles. 

7.1 Quantization of a field 

For the quantization of the fields, (j> r (x), the field and the conjugate fields, ir s (x), are promoted 
to field operators, i.e. 



These operators are Heisenberg operators and the following equal-time commutation relationship 
are imposed [34] : 




(7.1a) 




(7.1b) 




(7.2b) 



(7.2a) 



These commutation relations are imposed to ensure that the correct spin statistics are satisfied 
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Since the field operators are time-dependent their dynamics are determined by Heisenberg's 
equation of motion [35]. Heisenberg's equation of motion states that for an operator 0{t) the 
time-evolution of the operator is given by 

ih^ = [0(t),H], (7.3) 

where H is the Hamiltonian (16.6H . 

The field operators are defined in the space defined by the state vectors of the particles described 
by the field. To obtain an explicit expression of these operators they are Fourier decomposed in 
terms of a complete set of wave functions in the space [35] : 

4> r (t, x) = a ri (t) Uri(x) , (7.4a) 

i 

4>l(t,x) = Y,at i (t)u* ri (x). (7.4b) 

i 

The operator properties and the time-dependence of cp r (t,x) are carried by the expansion coef- 
ficients, a r i(t) [35] . 

7.2 Relativistic Scalar fields 

From Chapter [5j relativistic, spin-0 (scalar) particles are described by the Klein-Gordon 
equation. The Lagrangian density that describes a real scalar field, (f>(x), with particles of mass 
m, is (in natural units) [35j: 



L = m <p[x ) 

= h^{x)d^{x) - ^m 2 <j)(x) 2 . 

The equation of motion of the field is derived by applying the Euler-Lagrange equation (16. 4h 
and is 



d^4> + m 2 (j) = 0. 



(7.5) 
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The conjugate field to <f>(x) is denned analogous to the conjugate momentum of a classical field 
figa by 



7T = 



dC 



and the Hamiltonian density (|6.7[) is 



(7.6) 



= l(^ + (V0) 2 + mV). 

7.2.1 Quantization 

In canonical field quantization the following substitutions are made 

4>(x) — > 4>(x) , 

ir(x) — > tt(x) , 

and the operators have to satisfy the following equal-time commutation relationships: 



4>({t,x), TT(t,x') 

4>{t,x), 4>(t,x') 



iS(x — x) 

Tt(t,x), Tt(t,x') 







(7.7a) 
(7.7b) 



It is shown in Ref. that the scalar field has to satisfy the above commutation relations, 
otherwise the principle of microcausality will be violated. 



The Hamiltonian of the quantized field is simply found to be [35] 

H = ^Jdp(n 2 + (V4>) 2 + m 2 ^ 2 y (7.8) 

7.2.2 Expansion of operators 

The equation of motion of <fi(x) follows from Heisenberg's equation (|7.3p (keeping in mind 
that 7r(x) = <p(x)), 

4>(t,x) = (V 2 -m 2 )(j){t,x), (7.9) 

which is just the free Klein-Gordon equation (|5.2p . Expanding the operator in terms of the plane 
wave basis for the operator space, the wave function in Fourier decomposition of the operator 
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(I7.4j) becomes 

u p {t,x) = N p e ip x 

where N p is just a normalisation constant. Thus 4>{t,x) can be written as [35] : 

4>(t,x) = J dpN p e ipx d p (t). (7.10) 
Substituting the above equation back into Eq. (|7.9p . it is found, as shown in Ref. [35], that 

a P (t) = -(p 2 + m 2 ) d p (t) 
which means that solutions for d p (t) are given by 

a p (t) = 4 1J + a p 2) e +i ^\ (7.11) 

where oo p is 



oo p = \fp 2 + m 2 



and cip^ and a p 2 ^ are constants in time [35]. Since the original 4>{x) was a real- valued classical 
function, the operator should be hermitian = 0), and therefore one of the operators can be 
expressed in terms of the other as 



— w_ p 



Thus the free particle solution for the quantized scalar field (I7.10P becomes 

$(t, x) = J dp N p {a p e^-^vt) + a\ p e -i(p-*-"pt)^ ) (7.12) 

and the conjugate field operator, 7r = ((>, as 

7r(t, a;) = y dp N p (-ito p ) (a p e^-^vt) _ a t e -*(p-«"-«P*)^ . (7.13) 

As is shown in Ref. [35], by substituting the expressions for the field operator (|7.12p and the con- 
jugate field operator (|7.13p back into the commutation relations (|7.7p the commutation relations 
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for d v and d v can be obtained and are 



a. 




iS(p-p') 



(7.14a) 



a. 




= 0. 



(7.14b) 



Thus dp and a p satisfy the commutation relationship for creation and annihilation operators for 
scalar particles (bosons) [33J. Scalar particles are described by symmetric states. The states 
that the field operators will operate on are clearly symmetric from the above commutation 
relationships, since for an arbitrary state |3>) 



The Hamiltonian is constructed from Eq. (|7.8|) . by replacing the field operators by their 
expansions [Eqs (|7.12p and (j7.13j) ] . After some algebra (as is shown in Ref. [35] ) the Hamil- 
tonian is expressed as 



If the above Hamiltonian operates on any state, even the vacuum, it will have infinite energy, 
due to the contribution of the last term in Eq. (|7,15p . This contribution is attributed to the 
filled negative energy states that defines the vacuum [35]. To remove this constant contribution, 
Hamiltonian is defined to contain the normal-ordered product of the operators. Normal-ordering 
entails that all creation operators have to stand on the left and all the annihilation operators on 
the right [when operators are switched around the commutation relationships (|7.14p have to be 
taken into account] . The Hamiltonian, as the normal-ordered product of the field operators, is: 



oJ fit |$) = a£at |$) . 




(7.15) 
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7.3 Dirac fields 

To describe a field of Dirac (spin-^) particles the ansatz is that the Lagrangian density for 
the field is given by 

C = iijr {x)—ij){x) + ivpUx)ot ■ Vipix) — rmp^ (x)j3ip{x), 
ot 

where \{j(x) and the conjugate field ip^(x) [see (|7.17p ] are classical fields describing the Dirac 
particles with mass m. This Lagrangian density can be written covariantly using the definition 
of the 7-matrices (|5.1ip as well as the Dirac adjoint (|5.14p . as 

C = tft(x) [ij IM d^tp(x) — m] ip(x). 

Again the equation of motion can be derived using the Euler-Lagrange equation (|6.4p and is 

[i^df, -m] ip{x) = 0. (7.16) 

The conjugate field is 

fir 

tt = -j = itf (7.17) 

dip 

and the Hamiltonian density is 

n = n^-c 

= ip '(#)(— icx. • V + /3m)ip . 

7.3.1 Quantization 

The Dirac-fields are quantized by replacing the Dirac-spinors by field operators, 

ip(x) — > ip(x) , 
tjj\x) — > $ '(&), 

where the field operators satisfy the following equal-time anticommutation relationships, 

sc), $p{t,x')\ = 5 af3 5(x-x'), (7.18a) 
{$ a (t,x), M^ x ')} = {$(*,a:), $(t,sc')} =0. (7.18b) 
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In contrast to the quantization of the scalar field, the anticommutation relationships of the Dirac 
field operators are defined so that Fermi-Dirac statistics is obeyed [35] . 

7.3.2 Expansion of operators 

As is shown in Ref. [35J the quantized field operator satisfies the equation for the Dirac field 

i^(t, x) = {—la • V + /3m)$(t, x) , 

which is 

(i7 M d M -m)i>{x) = 0, (7.19) 
when expressed in a covariant manner. 



The operator for the Dirac field can be expanded in terms of the solutions to the free Dirac 
equation ()5.9|) . These solutions can be explicitly written (from Ref. [35]) as 

tf'ft") =(2^PVf W^(P,e "' < "''"' t> • <7 - 20) 



where r runs for 1 to 4. r = 1,2 denotes the positive energy solutions (E = uj p = \/ p 2 + m 2 ) 
and while the other two values denote the negative energy solutions (E = —oo p ). Thus e r = 1 for 
r = 1, 2 and e r = —1 for r = 3, 4. w r (p) are the Dirac spinors and from Ref. [35] they have the 
following orthogonality and completeness properties: 

wli(e r 'p) w r (e r p) = — 5 rr i , (7.21a) 

m 

w r '(p) w r (p) = e r 5 rr i , (7.21b) 

4 

^2w ra (e r p) wlp(e r p) = —Sa/3, (7.21c) 
r=l m 
4 

^e r io m (p) w r/3 (p) = 8 a p. (7.21d) 

r=l 
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The expansion of the field operator is thus given by 

4 

i>(t,x) = ^2 dp a(p,r) il$\t,x) 

r=l 

= ±Jj0m^ p ^^^>' < 7 - 22 > 

and that of the conjugate field operator by [35] 

ft(t,x) = Yl / dpa\p,r) 4 ] \t,x) 
r=l ^ 

It can be shown (as in Ref. [35]) that the operators a(p, r) and a\p, r) satisfy the same anticom- 
mutation relationships imposed on the field operators (|T. 18|) : 

\a(p,r), a}(p',r')} = 5(p-p')5 rr >, 
[a(p,r), a(p',r')} = {a f (p, r), a^p', r')} = 0. 

These anticommutation relationships describe anti-symmetric states which are occupied by fermions 
(spin-7^ particles) [33], since for an arbitrary state |<3?) 

a\q,r') a^(p,r) |3>) = -a\p,r) a)(q,r') |$) . 

The Hamiltonian of a quantized Dirac field is given by [35] 

H= j dp 1^2 w P a f (p,r) a(p,r) - J^w p a+^r) a(p,r) . (7.24) 

\r=l r=3 / 

This Hamiltonian can have an infinite (negative) value, for the last term adds the contribution 
of all the negative energy states (r = 3,4). In the Dirac interpretation of the negative energy 
states these states are all filled to form the vacuum. The only contribution of the negative energy 
states that should be considered is if an negative energy particle has been excited to a positive 
energy state, leaving a hole in the vacuum. To have a sensible interpretation of the Hamiltonian 
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the contribution of the vacuum has to be removed, i.e. 

. 2 4 

H = i dp\^u p at far) a(p,r) + ^uj p (l - a\p,r) afar)) 

J r=l r=3 

r 2 4 

= y dp^^Up a\p,r) d(p,r) + ^o; p (a fa r) a ] {p, : 



(7.25) 



r=l 



r=3 



For Eq. (|7.25p to be expressed as a product of normal-ordered operators, d(p,r) with (r = 1,2) 
are the annihilation operators for particles and a){p,r) with (r = 3,4) are the annihilation 
operators for antiparticles (holes). Similarly a^(p,r) with (r = 1,2) are the creation operators 
for particles and a(p, r) with (r = 3, 4) are the creation operators for antiparticles. To more 
easily distinguish between operators for particles and antiparticles, the spinors are renamed as: 



w x fa) 
w 2 fa) 

Wzfa) 
Wifa) 



u(p,+s) , 
ufa-s) , 
v(p, -s) , 
v(p, +s) , 



and the operators as: 



a fa 1) 
a(p,2) 
a fa 3) 
a fa A) 



b(p, +s 
bfa -s 
$(p,-s, , 
tf (p, +s 



where s denotes the spin projections. The expression for the field operators therefore changes to 
(from Ref. [35] ) 



and 



ft(t, x) = J 7^72 ^/J ( gt (Pi s ) "(p. e+i(p ' x) + rfA (p. s ) "(Pi s ) e " <(p,<B) ) • ( 7 - 27 ) 
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7.4 Vector boson fields 

Vector bosons are spin-1 particles and the commonly known ones are the massless photon and 
the massive and Z° bosons [35] . The u- and p-mesons are also spin-1 particles [21j . 

The quantization of the fields describing the massive uj- and p-mesons are of particular interest 
in this work, since these mesons are included in the description of nuclear matter studied in 
Chapters and EE The quantization of massive vector boson is to a large extent based on the 
quantization of the photon, a massless vector boson, and therefore the quantization of the photon 
field will be discussed first. 

As is explained by Greiner et al. in Ref. [35] the vector boson fields seem to be overdetermined 
since one would usually expect four degrees of freedom to be associated with a vector field (that 
of a scalar component, longitudinal and two transversal polarization directions). Photons, being 
massless, have only two observed degrees of freedom (two transversal polarization states), while 
massive vector bosons additionally also have a longitudinal polarization state. This reduction 
of the degrees of freedom has to be taken into account when quantizing the theory by imposing 
additional constraints on the field operators |35| . 

7.4.1 Photon field 

The electromagnetic interaction between particles via the exchange of photons is described 
by Maxwell's equations. (For a complete description of the electromagnetic interaction see for 
instance Introduction to Electrodynamics by D. J. Griffiths |15j.) Maxwell's equations can be 
expressed in a covariant manner by defining the field-strength tensor (or field tensor), F^ IU , 







E 1 



E 2 



E 3 



E 1 







B 3 B 2 



E 2 



B 3 







B 1 



E 3 



B 2 



B 1 







as 



(7.28a) 
(7.28b) 
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where j u = (p, j) is the electromagnetic four-current density. 

In terms of the electromagnetic vector potential, A^(x) = (A°(x), A(x)), the field tensor can 
be expressed as [55] 

= d»A v (x) - d v A»(x). 

Eq. (]7,28ap then leads to a second-order wave equation and field equation for A^ 

d^A v {x) - d u {d^{x)) = f{x), (7.29) 

from which it can be seen (by taking the four-divergence) that the electromagnetic current is 
conserved, i.e. d u j v (x) = 0. 

Since the addition of any scalar function A(t, x) to A^ in the form of 

A'^(x) = A»(x) + a^A(x) , (7.30) 

leaves F^ u invariant and thus also the field equation (|7.29p . Since Eq. (|7.30p describes a local 
gauge transformation, A^ is therefore gauge invariant |35j . 

The property of gauge invariance leads to technical complications in the theory |35| and an 
additional constraint ("gauge") needs to be imposed on the field to remedy the situation (in this 
case by reducing the degrees of freedom of the field) |34j . There are various different conditions 
that can be imposed on the photon field, one of which is the Lorenz condition, 

= 0. (7.31) 

By demanding that the Lorenz condition is satisfied, is known as the Lorenz gauge. (This con- 
dition is mostly misattributed to H. A. Lorentz [36].) 

The photon field equation (I7.29P can be derived from the following Lagrangian density [34] . 

C = - l -(d u A,WA»)-j,A» 
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and the conjugate field to the photon field is therefore 

dC 



d(d A») 



d u A L 



(7.32) 



7.4.2 Quantization 

The photon field is quantized by considering the field and the conjugate field as field operators 
and imposing the equal-time commutation relationships [35] 



A»(t,x),TT u (t,x') =%rf 6 3 (x-x'), 
A t *(t,x),A v (t,x')] = l^(t,x)X{t,x')} =0 



(7.33a) 
(7.33b) 



where , rf jV is the metric as defined in (12. 3p , 



Using the expression for the conjugate photon field (|7.32|) . Eq. (|7.33ap becomes 



A»(t,x),d°A v (t,x') 



irf w 5 3 {x-x' 



(7.34) 



The commutator of the divergence of the field operator, d^A 1 *, and the field operator can be 
evaluated, using Eqs (|7.33p and (|7.34p . and yields 



d^{t,x),A u {t,x') 



d A°(t, x) + V ■ A(t, x),A u (t, x') 



tr°{t,x),A u (t,x') 
ir] u0 5 3 (x-x') ^0. 



+ V 



A(t,x),A u (t,x') 



From the above it can be seen that d^A^ ^ and thus the canonical quantization procedure and 
the Lorenz gauge are not compatible |35j. 



If the Lorenz gauge is not imposed "scalar" photons can be created that arbitrarily add en- 
ergy to the system. Since only transversal photons are observed, another constraint has to be 
imposed to be able to quantize the photon field. In the Gupta-Bleuler method only the state 
vectors in the Hilbert space, |$), that satisfy 



= 0, 



(7.35) 
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are admitted. By imposing condition (|7.35p on the states in the Hilbert space the photon field 
can be quantized by demanding that the photon field operators satisfy the conditions given in 
CLMj) . [35] 

7.4.3 Massive vector bosons 

A massive neutral vector bosons (spin-1) field, A^(x), is described by the Lagrangian density 
[35], 

C = -\f^F^ + ^m 2 A fM (x)A^x)-j^x)A^x), (7.36) 
where F^ v is once again the field tensor, 

= d^A u (x) - d u A^(x) 

and jf_t(x) the current density. 

Using the Euler-Lagrange equation ()6.4p . the equation of motion of the field is calculated to 
be 

+ m 2 A v {x) - f{x) = 0, (7.37) 



which is the Proca equation |34| . 

By taking the four-divergence of the Proca equation, one finds that 

d v A v {x) = \d v f{x). 

If one assumes that there are no sources, i.e. j v {x) = 0, or that the source current is conserved, 
i.e. d v j v (x) = 0, the above condition on the field reduces to 

d v A v = 0, 

which is the Lorenz condition (|7.3ip . Thus the Proca equation for massive vector bosons auto- 
matically satisfies the Lorenz condition under the assumptions mentioned and simplifies to 

d^A u (x) + m 2 A v {x) = f(x). (7.38) 
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The conjugate field (17. 6p is derived from the Lagrangian density (17.36P as 

dC 



7T' 



(7.39) 



Since F^ v has no diagonal elements, 7r° would be zero and tt 1 = E % . Because the Proca equa- 
tion automatically satisfies the Lorenz condition the time-component of the conjugate field is 
a dependent variable. This is to be excepted since by automatically satisfying the Lorenz con- 
dition the Proca equation only has three degrees of freedom [35J. As a dependent variable the 
time-component of field is given in terms of the degrees of freedom (E l ) as 



A = — U • E. 



(7.40) 



7.4.4 Quantization 

The canonical quantization of the massive vector field is achieved by imposing the following 
commutation relationships: 



A\t,x),E\t,x') = -i5ij5 3 (x - x') 



A%x),A j (t,rf) = E i (t,x),E j (t,x') =0. 



(7.41a) 
(7.41b) 



It may seem that the quantization of the field might not be covariant, since it is only expressed 
in terms of the spatial components of the field operators. According to W. Greiner et al. a 
covariant quantization can be achieved by defining the time component of the field operator to 
be [35]: 

A = — • E. 

From (|7.41|) this operator satifies the following commutation relationships 



A(t,x),A° (t,x') 
A°(t,x'),A°(t,x') 



1 

m z 
0. 



V6 3 {x-x'), 



CHAPTER 8 
Quantum hadrodynamics 



8.1 Introduction 

As mentioned in Chapter H J. D, Walecka introduced a description of nuclei and nuclear 
matter based on the exchange of mesons in 1974. This description will be referred to as quantum 
hadrodynamics or QHD. 

Since nuclei and nuclear matter are complex systems, there exist various models of nuclei and 
nuclear matter of which QHD is one. All models need some form of experimental input to con- 
strain the model. In the case of QHD it is the coupling constants (or coupling strengths) between 
the different meson and nucleon fields that are the unknown parameters. These couplings are 
determined by fitting the calculated properties of nuclei and nuclear matter (such as the prop- 
erties that were described in Sec. I4.1j) to the experimentally observed (or inferred) values. In 
the last few decades various other parameter sets for QHD have been calculated by different 
authors by fitting various other observed properties of nuclei and/or nuclear matter. These dif- 
ferent parameter sets differ amongst one another since they might include other meson fields 
and/or different couplings between the various fields and/or recalculated values of the couplings 
constants. 

In Sec. 18.21 the basic philosophy and calculations will be introduced within the context of the 
original parameter set, QHD-I, as introduced by J. D. Walecka |24} I21j. 

In Chapter [9] other parameter sets and variations based on QHD-I will be discussed. 

8.2 Formalism 

Quantum hadrodynamics I (QHD-I), also known as the a - uj model, is the original and 
simplest parameter set of QHD. It models the nuclear force by the exchange of neutral (isoscalar) 
scalar sigma (a) mesons and neutral (isoscalar) vector omega (w) mesons. These mesons have 
been found to be the most important in describing the properties of nuclei and nuclear matter 
|25j . The baryons included in QHD-1 are protons and neutrons. The scalar meson gives rise to a 
strong attractive central force and a spin-orbit force in the nucleon-nucleon interaction, while the 
vector meson gives rise to a strong repulsive central force and a spin-orbit force (with the same 
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sign as the spin-orbit force of the scalar meson) |23j. Thus the nucleon-nucleon interaction is 
described by three fields, the baryon-, scalar meson- and vector meson-fields. No charged mesons 
are included in this parameter set (i.e. the electric properties of the baryons are not considered) 
and the masses of the proton and neutron are taken to be equal in QHD-I. 

The Lagrangian density of QHD-I (in natural units) is given by [21J: 



+ 2 



$(x)[^(id» - g v V»(xf) - (M - g s <t>{x))\i>{x) 
l -{d^{x)d^{x) - m 2 s <j> 2 {x)) - + l -ml V^x)V»(x) 



where 

• V denotes the vector meson field, 

• (ft denotes the scalar meson field, 

• m u and m s the different meson masses and M denotes the nucleon mass, 

• g v and g s are the vector and scalar coupling constants and 

• V^v = d^V^x) - dyV^x). 

The Lagrangian density is a Lorentz scalar and has the dimension of length . 

In constructing the Lagrangian density it was assumed that the neutral scalar mesons couple to 
the scalar density of the baryon field and that the neutral vector mesons couple to the conserved 
baryon current |21j . 



Using the Euler-Lagrange equation (|6.4p the equations of motion of the different fields can be 
derived: 

d ll d ll 4>(x) + m 2 s <f>(x) = g s ip(x)ip(x) (8.1a) 
+ mlV u (x) = gv^xW^x) (8.1b) 
l^-g v V^(x)) - [M - g s <t>{x))]i>{x) = 0. (8.1c) 



Eq. (|8.1ap is analogous to the equation of motion of the free scalar field, the Klein-Gordon equa- 
tion (|7.5p . with the baryon scalar density, vp(x)il>(x), as the source term. 
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Eq. (|8.1bj) is analogous to the Proca equation (I7.38P , with the coupling of the vector field to 
the conserved baryon current as the source term. 

The equation of motion of the baryon field is the Dirac equation (|7.16p with modified mass 
(due to the scalar field) and interaction with the vector field introduced in a minimal fashion 



Since the equations of motion (18. 1|) are non-linear, coupled equations they are very difficult 
to solve and are therefore approximated. The approximation that will be used is the relativistic 
mean-field approximation. 

8.3 Relativistic mean-field theory 

In other quantized theories, such as QCD, a perturbation expansion in the coupling strengths 
can be used to approximate complicated equations. Within the framework of QHD the cou- 
pling strengths are large and thus the higher-order terms in a perturbative approximation will 
diverge. Within the context of the study of dense nuclear systems an approximation that is 
increasingly valid as the density increases, is also preferred. In the relativistic mean-field approx- 
imation (RMF) within QHD, the meson field operators are replaced by their ground state (|3>)) 
expectation values, which are classical fields, in the following way |21j : 



The mean-field approximation (|8 . 2[) can be motivated if a system of B baryons (which are at rest 
and at zero temperature), occupying a box of volume V is considered. Since the system is static, 
the baryon flux, given by 4>(x)Y'4'(x), will be zero [22]. If the baryon density, B/V, is increased 
the source terms on the right-hand side of Eqs ()8. la|) and (|8.1bp will become large. If the source 
terms are large enough the meson field operators can be approximated by their ground state 
expectation values, instead of evaluating the operators at every spacetime point in the box. For 
a stationary, uniform system </>o and Vo will be constants, independent of space and time. The 
spatial components of (Vn) will vanish, since, as mentioned above, the baryon flux is zero and 
the system is at rest. Another illustration for the vanishing spatial components of (Vu) is given 







in Fief. 0). 
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In the mean-field approximation the meson field operators are replaced by their ground state 
expectation values in order to simply the solution to the field equations. However, the nucleon 
field operators remain operators. Since only the ground state expectation values of the meson 
field operators are considered in the approximation it follows that the baryon operators (the 
baryon sources to which the meson fields couple) must be evaluated by operating on the ground 
state. As such these baryon operators in the equations of motion of the meson fields (|8.ip are 
replaced by their normal-ordered ground state expectation values |37j : 



ip(x)ip(x) — >■ ($ \:ip(x)il)(x):\ $) = (ipip) 
'4){x)^il)(x) — > ($ \:%[){x)^%l)(x):\ 3>) = (VW) 



(8.3a) 
(8.3b) 



where the spatial components vanished due to the consideration of the ground state. The normal- 
ordered expectation value is taken since the contribution of the filled negative energy baryon 
states (the vacuum) is ignored, since the vacuum has a (infinite!) constant energy. Thus only 
the positive energy baryon states are considered. This is known as the no-sea approximation [42] . 



Since (j)Q and Vo are constants in the RMF approximation the equations of motion of the fields 
(1870) reduce to 



mlVo = g v (V>7oV'/ 
[i^-g v -f V -(M-g s (j)o)}^ = 0. 



(8.4a) 
(8.4b) 
(8.4c) 



The RMF QHD-I Lagrangian density is: 



C-RMF = i> 



il^ - g vl °V -{M-g. 



From the above Lagrangian density the RMF energy-momentum tensor of QHD-I can be con- 
structed using the definition of the energy-momentum tensor (|6.13p as well as keeping Eq. (|8. lcj) 
in mind, and is given by 

(T^)rmf = iH 1 ^ - VT ( " \ ml4 + \ mlV 2 ) . 
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The energy density and the pressure, using (16. 14|) for QHD-1 is |21j 




{iipj d ip 
(iipjodoip) 



( - 2 m ^o + 2 m - F o 

1 2/2 1 2 T /2 

2 m s ( t , o ~2 m u V o ' 



2 



) 



+ 



(8.5a) 



P 




(8.5b) 



8.4 Evaluation of expectation values 

Since ip and ip (the baryon field operators) are still operators, their expectation values have 
to be calculated to obtain sensible information for, amongst others, the energy density and the 
pressure. As mentioned previously, it is the ground state expectation values of the baryon field 
operators that have to be evaluated. 

These expectation values can be evaluated by constructing an explicit expression for ip. For 
the RMF approximation a uniform static system is assumed, so ip is expanded in terms of the 
momentum eigenstates [21]. The momentum eigenstates are denoted by k, which is the wave 
vector. The wave vector is related to the momentum (p) through the expression of the de Broglie 
wavelength [38J, 




In natural units the momentum and the wave vector are equivalent. Solutions for the field 
operator are thus of the form [analogous to the free single-particle solution for a Dirac particle 



where ip(k, s) is the four component Dirac spinor (s denotes the spin index) and e{k) the energy 
associated with specific state denote by k [21]. Substituting Eq. (|8.6f) into Eq. (j8.4cj> yields 




p = hk . 



ip(x) = ip(k, s) e 



: ik-x—ie(k)t 



(8.6) 



( - Jitf + 7o e(fc) - g vlQ Vo - (AT - 9s ^ ))iP{k, s) = . 
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By multiplying the above equation by 70, it reduces to (having reverted back to the notation of 
the a- and /3-matrices) 

(p 2 (e(k)-g v V )^(k,s) = (a-k + Pm*^(k,s), (8.7) 

where 

• m* is the reduced nucleon mass, defined as 

m* = (M — g s ^) . (8.8) 

The presence of the scalar meson field thus shifts the mass of the nucleons |21| . 

Squaring both sides and keeping the properties of the a- and /3-matrices (|5.6p in mind one 
finds that 

e ± (k) = g v V ± Vk 2 + m* 2 . (8.9) 

From Eq. (|8,9p is it clear that positive and negative energy solutions, which are symmetric 
around g v Vo, are still found. Defining positive and negative energy spinors u(k,s) and v(k,s) 
that in general satisfy the same orthogonality and completeness properties as in (|7.2ip . as well 
as the following normalisation condition [21] 

«^(fc, s)u(k, s') = v\k, s)v(k, s') = 5 ss i , 

the general solution for the baryon field operator (|8.4cp can be written as [21] 

=E / (^7^( S(fc ' s) n(fe ' s) e lk - x - ie(+Hh)t + S(k,s) v(k,s) e-ik-'-iJ-Hk)^ 7 (8 .io) 

where the particle and anti-particle creation and annihilation have been included (using the same 
notation as in Sec. I7.3p . 

The ground state |<3?) is defined by kp, called the Fermi-wavenumber: kp denotes the filled 
positive energy baryon momentum state below which all positive energy states are also filled. 
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The ground state in the the no-sea approximation has the following properties: 

d(k,s) |$) = 

St(fc,«)|$) = 

b(k,s) |$) = 

The expectation value in the expression for the RMF value of Vq (|8.4bp can be evaluated by using 
the construction of ip{x) f|8. 10|) and taking the normal-ordered expectation value, as in Eq. (|8.3j) . 
In normal-ordering the expression the following should be kept in mind: The anti-commutation 
relationship of b(k, s) and dJ(k, s), which is the same as the creation and annihilation operators 
of the general Dirac field given in Sec. 17,31 as wen as the- properties of the ground state. Then 
it can be shown that, after some algebra, (ifr'tp') is given by 




(8.11) 



where 7 is the nucleon spin-degeneracy of the state. From (|8.11|) it is clear that ("^'V) 1S the 
nucleon (baryon) number density at |fc| = kp, which will be denoted by p. 

From (|8.5|) the expressions for the energy density and the pressure for QHD-1 in the RMF 
approximation become 

e = ^(-ia-V + f3m*+g v V )^-^mlV 2 + ^m 2 J 2 , (8.12a) 
P = l^(- ia - V)v)-^m2y 2 + im^. (8.12b) 

The expectation values in (I8.12P can be evaluated in the same way as the one in Eq. (18. 11 f) . 
After some algebra [using in particular the properties of the spinors (|8.7p as well as Eq. (|8.9p ]. 



V \k\ 

V |fe| < kp 

V Ifcl > k F 
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the expectation values in the expressions for the energy density and the pressure reduce to: 



The remaining expectation value, that of vipip) in the equation of motion of the scalar meson 




field (|8.4ap . can be evaluated in this explicit manner, but was not done in this work. Rather 
other considerations (as explained in Ref. [21]) will be employed to yield the expectation value 
(see below and Sec. I8.6|> . 

Another method to evaluate the expectation values is given by N. K. Glendenning in Ref. [3]. It 
states that the expectation value of an operator (r) in the ground state can be given in terms of 
the expectation value of the single-particle state, 





k.s ' 



where 



• k denotes the momentum, and 



• s the spin of the single-particle state. 



The expectation value in the many-nucleon system is then [2] 





where 



• e(k) = g v Vo + (k 2 + m* 2 )^ is the positive single-particle energies (18. 9j) . since only the 
ground state is considered, 



• fj, is the chemical potential/Fermi energy (single-particle energy e(kp) ), and 



• Q[fi — e(k)] is a step function with 



G[/i-e(fe)] = < 



1 if \k\ < k F 



if |fc| > k F 
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All the possibilities for T will in general appear in the Dirac Hamiltonian [4]. The Dirac Hamil- 
tonian Ho can be constructed from the equation of motion of the nucleon field (18.4ch and the 
expression of ij) (18, 6h as 



H D = ■ k + g vl0 V + m*] . (8.15) 
Taking the single-particle expectation value of Hp 

= e(k) 

k,s \ J k,s 

= e(k), (8.16) 

where (ify 'if) fc = 1 (this can be seen from the expression for tp and the normalisation of the 
spinors). Taking the derivative of the left-hand side of Eq. (|8.16|) with respect to any variable, 
C, yields 

as ip(x) is an eigenfunction of Hp, the second term on the right vanishes [4]. Thus Eq. ([8.170 
[considering Eq. f|8.16|> ] yields 

Any expectation value can therefore be obtained by using the general expression ([8.140 : taking 
the derivative of Hp for the appropriate choice of £ to yield T and obtaining the expression for 
the single particle expectation value of T from Eq. ([8.180 . 



Using this method the same results are obtained as with the first method, just with much less 
effort. The expectation value in the expression for the energy density ([8.12ah . when V operated 
on ifi ([Q|) . using Eqs flBZ!> and (j8l)]) . is: 

(V (a • k + /3m* + g v V )ipj = (V k 2 + m* 2 + g v V^j ip\ 
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and using the value of the single-particle expectation value of Hjj (I8.16p . the expectation value 
in the energy density yields 

= E / J^(y/k 2 + m^ + g v V )e[fi-e(k)] 

= 9vV p+—^ [ dky/W 
(27rr Jo 

The expectation value in the expression for the pressure (|8.12b|) is 



m* 2 . 



= (V>(7 ■ k)ip) , 
and the single-particle expectation value of the operator 7 • k is 

J 9Hd i\ u 

k using ^M) 



dk 
kk 



Vk?~ 



m 



71 ' 



and using the general expression for the expectation value (|8.14[) 



T /"^ fc 2 

/ rffc =. (8.19) 



(2vr) 3 7 Vfc 2 + m*' 2 
A result for (V^) can a l so be obtained in this manner: taking M (the nucleon mass) as £, the 
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-^proton -^neutron m s m ui 
939 939 520 783~ 

Table 8.1: Particle masses (in MeV) for QHD-1 used in this work. 

single-particle expectation value is 

/ , f J.) 

^ dM * k > s 

d 

e{k) using (f87T8^ 



dM 
m 



Vk 2 + m* 2 



Using the general expression for the expectation value (|8.14p and the above result, (ipip) yields 

= \ r- ~, (8-20) 

2tt 2 Jo ^k 2 +m* 2 

8.5 QHD-I parameter set 

The values of the particle masses used in this work for QHD-1 was taken from Ref. [22\ and 
are listed in Table [H7T1 



The values for the coupling constants used in this work for QHD-1 were taken from Ref. [22] 
(these values are preferred to those given in Ref. j21| . as those in Ref. [22] give a better 
description of the saturation properties of nuclear matter). The couplings constants are given in 
terms of C s and C v , which are defined as 

In Table 1531 the values for C s and C v and the corresponding values for g s and g v , using the values 
of the particle masses in Table 18.11 are given. 





r 2 


9i 


9v 


357.4 


273.8 


109.6 


190.4 



Table 8.2: Coupling constants of the QHD-I parameter set. 
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8.6 Equation of state 

Using the results for the expectation values obtained in Sec. 18,41 the RMF equations of motion 
for the fields f|8.4j) reduce to 

^o = 4 W = —2T~J kF dk 7^= 2 ■ ( 8 - 21a ) 
A. A/,U\ - _^ 



^ = -^rW) = ^P, (8.21b) 

while the energy density and pressure, from (|8.5[) and using Eq. (I8.21bjl to simplify the expression 
for the energy density, are 

11 'Y 

e = 2 m ^0 + 2™Zy« + pjjs / <* V* 2 + (M - s,*,) 2 , (8.22a) 
P = -imM + ^tf + \ (p^ f ^^=|=_) . (8.22b) 

The expression for <f>Q could have also been obtained from the thermodynamic argument that a 
closed, isolated system such as the one assumed for the RMF approximation will minimise its 
energy. Therefore an expression for m* (and therefore 4>q) can be obtained by minimising the 
energy density (|8.22ap with regards to m*. In this way the expression for m* is obtained, 

m * = M-^-^ f F dk^ =JlL= , (8.23) 
which is equivalent to Eq. (|8.21ap when the definition for m* (|8.8p is considered. 



Imposing spherical symmetry, the energy density (|8.22ap and the pressure (|8.22bp can be re- 
written as: 

e = -m^ + -raX + A / dk k 2 Vk 2 + m* 2 (8.24a) 
2 2 2tt z J 

P = --m 2 J 2 + - m lV 2 + I ( J- f kf dk-=¥=) . (8.24b) 

8.7 Observables of nuclear matter in QHD-I 
8.7.1 Symmetry energy 

To calculate the symmetry coefficient, £14, given by (|4.ip 
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the energy density must be expressed in terms of t. This can be done by considering the 
expressions for the proton (p p ), neutron (p n ) and baryon (p) number densities. These are given 
in terms of the different fermi momenta (k p for the protons and k n for the neutrons) of the species 
as 

1 f kp „ , 2 _ k p 



p p = — dkk 



IT 2 Jq 3lT 2 ' 

i r kn . 2 ki 



Pn = —r\ dkk 



IT 2 J 3ir 2 ' 

Pp + Pn- 



In the case of symmetric nuclear matter (t = 0) the proton and neutron fermi momenta are 
equal, k p = k n = Uf, and the baryon density reduces to 

p = 2x-Ly o dkk* = 2 ^. (8.25) 

The only terms in the energy density (|8.24ap that will contribute to 04 are the integrals over the 
filled baryon states. By defining e s as the isospin dependent part of the energy density, which is 



1 f 1 f p /— 

€fs= dkk 2 Vk 2 + m* 2 + — / dkk 2 Vk 

7T 2 JO 7T 2 JO 



2 _|_ m *2 ^ 



where the nucleon degeneracy, 7, has been explicitly included (as 2 for each baryon), the 
expression for 04 ()4. 1[) reduces to 

r^v 

.26) 



1 d /ll f " 1 1 f p \ 

- tt^t 0/ dkk 2 Vk 2 + m* 2 H ^ dk k 2 V k 2 + m* 2 

2 dt 2 \ P tt 2 J pn 2 Jo J 



t=o 



To evaluate the expression for 04 the proton and neutron fermi momenta are expressed in terms 
of t, i.e. 
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and 



= |^l^(l-t)V . (8.28) 



After substituting in the expressions for the different fermi momenta [ (|8.27j) and (|8.28|) ] into the 
expression for 04, as well as some algebra and remembering that k n and k p are equivalent at the 
evaluation point, (|8.26p yields 



a 1 



6\/k F + m* 2 



8.8 Summary 

The aim of this chapter was to introduce the basic concepts of QHD and RMF within the 
context of the QHD-I parameter set. This chapter will serve as the foundation for the derivation 
of the equation of state of nuclear and neutron star matter in other, more advanced QHD 
parameter sets. These parameter sets will be discussed in the next chapter. 



CHAPTER 9 
Advanced QHD parameter sets 

9.1 Introduction 

The coupling constants of QHD-I were fitted to reproduce the saturation of nuclear matter 
at a fermi wavenumber of 1.30 fm _1 and binding energy of -15.75 MeV [22]. Unfortunately this 
parameterisation produces too high a value of the nuclear matter compressibility (K). In 1977 
J. Boguta and A. R. Bodmer proposed that self-couplings of the scalar meson field should be 
included in the Lagrangian density and the coupling constants re-adjusted to reproduce K more 
accurately [39] . These self-couplings of the scalar meson field are included in all other QHD 
parameter sets studied in this work. 

Another addition is the inclusion of the charged (isovector) vector rho meson triplet (p°, p^) in 
all parameter sets discussed below. Since protons and neutrons practically only differ in terms of 
their isospin projections, the rho mesons are included to distinguish between these baryons and 
to give a better account of the symmetry energy [4]. (As these vector meson are charged, the 
reaction between a rho meson and a proton will differ from the reaction between a rho meson 
and a neutron.) 

The rho mesons are introduced in a similar fashion as the other mesons in QHD: the free La- 
grangian density describing the vector field is included in the QHD Lagrangian density as well 
as a coupling between the rho meson and the conserved isospin density [H [21] . 

The NL-SH parameter set is "an expansion of the ideas of Boguta and Bodmer" |40j and 
was formulated to reproduce the neutron radii of neutron-rich nuclei more accurately. This was 
to be achieved by fitting the coupling constants to reproduce the experimentally observed values 
of the root-mean-square neutron radii of neutron-rich nuclei [40 J . 

In 1994 Y. Sugahara and H. Toki introduced a self-coupling in the omega meson field to the 
Lagrangian density and formulated the TM1 and TM2 parameter sets. This was done to obtain 
better agreement with other models for nuclei and nuclear matter. The parameters were fitted 
to obtain agreement with the observed properties of unstable nuclei [41j . The TM1 parameter 
set is obtained when the fitting procedure is done to obtain agreement with light unstable nuclei 
and the TM2 parameter set for heavy unstable nuclei. 
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According to Ref. [32] the NL3 parameter set was formulated by G. A. Lalazissis et al. to 
(amongst others) further improve the predicted value of the compressibility, K. The NL3 
parameter set is based on the QHD-I Lagrangian density with the inclusion of rho mesons 
and self-couplings of the scalar field. This parameter set aimed to achieve a better descrip- 
tion of nuclei with large isospin asymmetry, i.e. nuclei far from the valley of beta stability. The 
parameters were calculated by fitting the predicted values of various properties (such as binding 
energies, charge radii, and neutron radii of spherical nuclei) to the observed values of nuclei that 
has a large isospin asymmetry [42] . 

As a further refinement to the TM1 parameter set, the PK1 parameter set was obtained by fitting 
properties (such as binding energies and charge radii) of a wide range of heavy nuclei. This was 
done to, amongst others, better reproduce the nuclear symmetry energy and compressibility [43] . 

B. Todd-Rutel and J. Piekarewicz introduced a coupling between the omega meson field and 
the rho meson field in an attempt to better describe the density dependence of the nuclear sym- 
metry energy without changing the saturation properties of nuclear matter |44| . The FSUGold 
parameter set was obtained by fitting the binding energies and charge radii of a variety of magic 
nuclei to the calculated properties. 

9.1.1 Photon field 

In general the photon (Coulomb) field is included in all the models mentioned above when 
studying the properties of finite nuclei. In this work the properties of infinite nuclear matter, 
specifically those which are applicable to the study of neutron star matter, will be studied. Since 
the Coulomb forces are much stronger than the gravitational forces on the nuclear scale, neutron 
stars are assumed to be charge neutral [3], and therefore the contribution of the photon field is 
ignored in this treatment of nuclear matter. 
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9.2 Lagrangian density 

The most general Lagrangian density, that describes the QHD-I, NL-SH, NL3, FSUGold, 
PK1, TM1 and TM2 parameter sets in QHD, is [37]: 

C = ^(x)[ 7 M H-<7*W) -Vf-r-b^x)- (M - g s cf>(x))]ijj(x) 

+ ±d»m(x) - \mtf(x) - |M(*)) 3 - ^M(*)) 4 



4! 



+ \rnlV»{x)V,{x) + £ {glV»{x)V,(x)) 2 



+ a„( 5 ^(*)v^))(^(x) . Mx) ) 



(9.1) 



where the field tensors have been defined as 



iiv = d^V v {x) - d u V^(x) 



dj) v (x) - d v bJx) 



with 



ip(x) the isodoublet baryon field (only protons and neutrons are considered in this work), 
4>(x) the sigma (scalar) meson field, 
V^(x) the omega (vector) meson field, 

b^(x) is the Lorentz vector field denoting the three isospin components of the rho meson 
fields, 

b" = (6?,6M. 

The charged rho meson fields (p^) can be constructed in terms of the first two components 
of W as @] 



1 



(9.3) 



and, 



r = (ti,T2,T3) the isospin operator. This operator is described in terms of the Pauli 2x2 
spin-matrices (|5.8f) as 

<T 

(J 
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It should be noted that since the nucleon field, if) is an isodoublet and therefore 



where is the tp p is the proton field and ip n is the neutron field, each consisting of the 4x1 
Dirac spinors, that the total dimension of ip is 8x1. Therefore r is in actual fact given by 



1 4 
U 



(9.4) 



and therefore the explicit expression for T3, using explicit representation of 03 (|5.8c|) . is 



1 




1 4 




1 4 


-1 




1 4 




-1 4 



(9.5) 



Using the Euler-Lagrange equation (16, 4h the equations of motion of the different fields are: 

A 



d^cj, + mi(f> + -gffl + 

.2 i v + C_ 

3! 



4 a3 



2 p S r ■ 3 |. 



d^ v + my + 2h v g 2 v V»V„g 2 .b v 



9p u 
—t ■ b„ 

2 ^ 



(M-g s 



g s ipip 

if, 
2 

0. 



(9.6a) 

(9.6b) 
(9.6c) 
(9.6d) 



9.3 Conventions 

In the literature two different conventions for expressing the Lagrangian density of QHD are 
used. These differences are obvious when the Lagrangian density and the values of the coupling 
constants of the NL3 parameter set in Ref. [H] are compared to those in Ref. 



The way in which the Lagrangian density ()9. 1 j) is expressed will be referred to as the Walecka 
convention, which is used in the original formulation of QHD by J. D. Walecka (see Ref. [21J). 
The other convention of expressing the Lagrangian density is mainly used by authors from Eu- 
rope and the Orient. This formulation will be referred to as the Ring convention. 



2 see IP 
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The differences are: 

1. In the Walecka convention the scalar field 4> is explicitly positive, where as in the Ring 
convention it is negative. 

2. The first four components in the nucleon field, tp, of the Walecka convention refer to the 
proton field and the last four to neutron field, whereas in the Ring convention the first four 
components of the nucleon field refer to the neutron field and the last four to the proton 
field. 

3. In the Walecka convention the coupling between the rho-field, b p and the isospin density, 
g p has an explicit factor of a half and in the Ring convention this factor is implied. 



Property 



Walecka 



Ring 



Reduced nucleon mass (m*) 

Coupling between the rho field and 
nucleon fields, g p 

Third-order self-coupling constant 
in the scalar meson field 

Fourth-order self-coupling constant 
in the scalar meson field 

Self-coupling constant in the omega 
meson field 

Coupling between the rho- and the omega 
meson fields 



Nucleon field, tp 



m* = M - g s 4> m* = M + g s <fi 

i)^(g p /2)T ■ b p ip ^rfg p T ■ b p ip 

« 92 
A #3 
C c 3 



A, 



Not applicable 

/ r 

1p = 



Table 9.1: Differences in Walecka and Ring expressions of the Lagrangian density 
of QHD parameter sets for nuclear matter. 

In this work the Walecka convention will be used. The expressions in the Ring convention 
will be given for certain properties, but will be identified as such. Both conventions are included 
in this work, because the different QHD parameter sets were studied in the specific context that 
they were defined in. It should also be noted that while certain terms will be present in one 
convention it might not be present in the other. This is due to the fact that different parameter 
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sets include different couplings and as such certain couplings might only be used in parameter 
sets studied by authors using a specific convention. 

It should also be noted that the coupling constants are defined separately for each conven- 
tion. 

The differences between the two conventions are contrasted in Table 19.11 

9.3.1 Note on m* 

The reduced mass (|8.8p . 

m* = (M — g s <j> ) , (9.7) 

is only defined in terms of the generic baryon mass, M. The baryons considered in this work are 
protons and neutrons and their masses are usually taken to be equal (939 MeV). However, the 
PK1 parameter set distinguishes between the proton and the neutron mass (see Table [93]) . This 
distinction should be taken (and is taken into account in this work) when calculating m* in the 
PK1 parameter set, but it is not explicitly shown in the expressions for m*. 

9.4 Relativistic mean-field theory 

As in Chapter [8] the equations of motion for the different fields are solved in the relativistic 
mean-field approximation, which assumes that the ground state of a uniform, static, spherical 
symmetric system is described and therefore the meson fields can be replaced their expectation 
values [37] : 

<p{x) > (cf>(x)) = (j)Q 

V»( x ) — > {V t (x))=g^ D V 
W{x) — ► (bZ(x))=g»°5 a3 b 

The spatial components of b^ and vanish due to the consideration of a uniform, static ground 
state (as was explained in Chapter |5J) . The ground state is also assumed to have definite spin 
and parity. Since the first two components of b^ can be written in terms of raising and lowering 
operators of the charged rho meson fields (|9.3j) . only the third component (that describes the 
neutral rho meson, p°) has a non- vanishing expectation value in the RMF approximation [3]. 
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Correspondingly the baryon sources in the equations of motion are replaced by their normal- 
order expectation values in the RMF grondstate, |<&), 



, tp(x)^(x) 



(<3?| :ij){x)i>{x): |$) = (x 
($| -4{x)^ip{x): |$) = (^7>) 



In the RMF approximation, the equations of motion of the fields (|9.6p reduce to 



9s<t>o 
9vV 

g P bo 




9v 
ml 

mA 



V-V) - UdvVof - 2A v (g v V )(g p b f 



V>W) -2A v (g v V ) 2 (g p b ) 



In the Ring convention^ the equations of motion are given by: 



(9.8a) 
(9.8b) 

(9.8c) 
(9.8d) 



m^ 



-5s (i>i>) ~ g2<t>l - 9300 

g P (^^3^) 

[i^d p - g vl °V Q - g p T 3 b - (M + g s 



(9.9a) 
(9.9b) 
(9.9c) 
(9.9d) 



In the RMF approximation the Lagrangian density (|9.ip reduces to 



Crmf = ^{id p - g v V (x)) - ^T 3 b — (M — 5^0) ^0*0 



1 

2 
1 

2 p 



4! 



+ T^p^o + A f {gvV ) 2 (g p b ) 2 



(9.10) 



3 The coupling between the rho- and the omega meson fields, A„, is not included in any of the parameter sets 
studied within the Ring convention. 
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Using the RMF Lagrangian density and (|6.14l) [keeping Eq. (|9.8d|) in mind] the energy density 
and the pressure are given by 



<T 00 > 

(i^odoip) - (C) , (9.11a) 



P = \{T U ) 



\ (iWdtf) + <£> , (9.Hb) 



where (£) is given by 



+ ±mlb 2 + A v {g v V ) 2 {g p bo) 2 . (9.12) 



9.5 Evaluation of expectation values 

The expectation values in the equations of motion of both conventions [Eqs (j9.8|) and (|9.9p ] 
and the equation of state (I9.1ip . were evaluated using the method given by N.K. Glendenning 
in Ref. [1], which is explained in Sec. 18. 41 



The calculation is essentially the same as in Sec. l8.4( with the only difference being that the 
expression for the Dirac Hamiltonian, Hp, and the positive single-particle energies (e(fe)J are 
slightly different. In this case, Hjj is given by 



Hd = 7o 



7 • k + g v j Vo + ^g P t 3 b + m* ( 9.1:]) 



and e(k) by 



e(k) = g v V + ^g P r 3 b + Vk 2 + m* 2 (9.14) 



where m* = (M — g s <fto) and the expansion for ip is taken to be as in (j8.6fl . with only the value 
for e(k) changing. 



Using the above it can be established that the expression of ("0 V ! ) wm n °t change from that 
in Chapter El except in this case the contribution of the baryons are evaluated separately (each 
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baryon has 2 spin states, and k n and k p denote the fermi momenta of the different species), 



E 



(2tt) 3 



(2tt) 2 



— * 3 + 

Pp + Pn- 



3tt 2 " 



dk 



(9.15) 



The same applies to the expression for (ipip), 



7T 



dk 



k 2 m* 



+ 



rrr 



-n-2 , 
7T ./0 



dk 



k m* 



(9.16) 



mr 



To evaluate OipW 3 ip\, a slightly different approach will be used. Since the baryon field can be 
expanded in terms of the different baryon components (in this case protons and neutrons), each 
described by a four component Dirac spinor, ip can be written as 



T3 is given by Eq. ()9.5p and therefore the effect of the operation of T3 will be to separate the 
contribution of the baryons into a positive (proton) and negative (neutron) parts. The single 
particle expectation value will therefore be given, since the states are normalised, by 

(y,t T3 ^) = { ^p T ^p\ = 1 

{ipt,T 3 ip n ) n = -1 



and thus 



_Lz~ 3 L_t 3 

Pp ~ Pn- 



(9.17) 



Since the following expression still holds 

de(k) 
dk 



■ k 



k ■ k 
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the expectation value in the expression of the pressure (I9.11b|) does not change, and is therefore 
still given by (IHTTOI) 

1 l' kp k 2 1 f kn k 2 

(iWdn/j) = dk - = + — / dk — =. (9.18) 

However, the expression for e(k) differs from the one in Chapter [8] [compare (|8.9p and (|9.14l) ]. 
The expectation value in the expression for the energy density is now given by 

<#7o^> = T^J J^e(k)e[p-e(k)} 

/dk 1 
7^3 {9vV + ~g p r 3 6 + V k 2 + m* 2 ) 6[/i - e(fc)] 

= ft>Vbp+ ^gpbo(pp-pn) +Y1 J j^p^vk 2 + m* 2 Q[p-e(k)} 

= gvV p + ^g p b [p p - p n ) 

1 f kp i 1 f k " i 

+ — 5 / dky/k 2 + m* 2 + — 5 / dkVk 2 + m* 2 (9.19) 

4tt 2 y 4vr 2 y 

9.6 Parameter sets 

The parameters for the QHD-I calculations done in this work are those given by B.D. Serot 
and J.D. Walecka in Ref. [22J. 

The Ring NL3 parameters are taken from the original article by G.A. Lalazissis et al. [32], while 
the parameters in the Walecka convention are from Ref. [44J. 

The parameters for the PK1 model are taken from the article by W. Long et al. |43| . while the 
NL-SH and TM1 (as well as TM2) paramater sets are taken from the article by M.M. Sharma 
and M.A. Nagarajan [40J and the article by Y. Sugahara and H. Toki [4T|, respectively. 
The parameters for the FSUGold model and the NL3 parameters in the Walecka convention are 
given by B. Todd-Rutel and J. Piekarewicz in Ref. 



The values for the coupling constants and masses for models studied in the Walecka conven- 
tion are given in Tables 19.21 and 19.31 respectively, while those for the Ring convention are given 
in Tables E31 and [931 
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Model 


g'i 


9v 


4 


K (MeV) 


A 


c 


At, 


NL3 


104.3871 


165.5854 


79.6000 


3.8599 


-0.01591 


0.00 


0.00 


S271 


81.1071 


116.7655 


85.4357 


6.6834 


-0.01580 


0.00 


0.00 


Z271 


49.4401 


70.6689 


90.2110 


6.1696 


+0.15634 


0.06 


0.00 


FSUGold 


112.1996 


204.5469 


138.4701 


1.4203 


+0.0238 


0.0600 


0.0300 



Table 9.2: Coupling constants of different parameter sets denned in the Walecka 
convention. All coupling constants are dimensionless, except for k which is given in 
MeV. 



Model 


-^proton 


-^neutron 


m s 


m w 


m p 


NL3 


939 


939 


508.1940 


782.5 


763 


S271 


939 


939 


505.000 


783 


763 


Z271 


939 


939 


465.000 


783 


763 


FSUGold 


939 


939 


491.5000 


783 


763 



Table 9.3: Particle masses (in MeV) of the parameter sets defined in the Walecka 
convention. 

9.7 Equation of state 

The equations of motion of the different meson fields ()9.8[) reduce to the following expressions 
when the results of Sec. 19.51 are taken into account. In the Walecka convention the equations of 
motion f|9.8j) reduce to 




(9.20a) 
(9.20b) 



g P b = %U(pp-Pn)-2A v (g v V ) 2 (g p bo)) (9.20c) 

Tibs* « £ ' 



Model 


g s 


9v 


g P 


52 (&Q v ) 


53 


C3 


NL3 


10.217 


12.868 


4.474 


-10.431 


-28.885 


0.00 


PK1 


10.3222 


13.0131 


4.5297 


-8.1688 


-9.9976 


55.636 


TM1 


10.0289 


12.6139 


4.6322 


-7.2325 


0.6183 


71.3075 


TM2 


11.4694 


14.6377 


4.6783 


-4.4440 


4.6076 


84.5318 


NL-SH 


10.444 


12.945 


4.383 


-6.9099 


-15.8337 


0.00 



Table 9.4: Coupling constants of different parameter sets defined in the Ring con- 
vention. All coupling constants are dimensionless, except for g<i which is given in 
fm . 
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Model 


-^proton 


-^neutron 


m s 


m w 


m p 


NL3 


939 


939 


508.194 


782.501 


763.0000 


PK1 


939.5731 


938.2796 


514.0891 


784.254 


763 


TM1 


938 


938 


511.198 


783.0 


770.0 


TM2 


938 


938 


526.443 


783.0 


770.0 


NL-SH 


939 


939 


526.059 


783.00 


763.00 



Table 9.5: Particle masses (in MeV) of the parameter sets denned in the Ring 
convention. 

The equations of motion of the fields in the Ring convention (]9.9[) are, after also considering the 
results of Sec. 19.51 as well as the difference between the conventions (Table I9TT]) . 



™%<f>0 = -M~2 / dk /,o = , = ? + ~2 / dk ~ 

V7r Jo Vk 2 + m* 2 n Jo Vk 2 + m* 2 ' 
-924-934 (9.21a) 

mlV = g v ( p n + P P ) - c 3 V 3 (9.21b) 

m l b o = 9 P {Pn-p P )- (9.21c) 

Using the expression for (£) ()9. 12[) and the results of Sec. 19.51 the energy density ()9.11ap and the 
pressure (|9.11bp in the Walecka convention yields 

- 2 m p b o ~ ^v(9vV ) 2 {g p b ) 2 + g v Vo{p n + P P ) + ^g P bo(Pp - Pn) 

1 f' kp 1 f' kn 

+^ / dk k 2 \Jk 2 + m* 2 H — = / dkk 2 Vk 2 + m* 2 (9.22a) 
Jo vr 2 J 

P = -\ml4- |(^ ) 3 - ^of + \m 2 X + ^{g v V f 

+\rn 2 p b 2 + A v (g v V ) 2 (g p b ) 2 

1 r kp k 4 1 f kn k 4 

H k / dk — = + — 5 / dfc . =■ (9.22b) 

3?r 7o V A; 2 + m* 2 37T Z Jo VfcS + m *2 



£; 2 m* 1 f kn „ /c 2 m* 
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In the Ring convention the energy density and the pressure are given by 

1 2/2 , 92 ,3 93 ,4 1 2 TA 2 c 3 T ,4 

e = 2 m ^° + T ^° + ~ 2 " ~ 7 

1 ,2.2 



-2 m p b + 3«Mj(Pn + Pp) + gpbo(Pn ~ Pp) 

1 f kp 1 i' kn 

H — ^/ dkk 2 \Jk 2 + to* 2 + / dkk 2 Vk 2 +m* 2 (9.23a) 



1 2,2 92 ,3 93 ,4 . * 2 T/ 2 , c 3 T ,4 , 1 



1 f fcp A: 4 1 /* fcri k 4 

+ — =■ / dk - + — =■ / dk - . (9.23b) 

3vr 2 7 VP+m* 2 3tt 2 7 Vfc 2 + to* 2 V 7 

9.8 Nuclear matter observables 
9.8.1 Symmetry energy 

The symmetry energy coefficient for the different parameter sets is calculated in exactly the 
same way as it was done for QHD-I in Sec. l8.7.iT except that in this case the energy density has 
more isospin dependent terms. 



The general expression for the energy density is given by equation (|9.22ap . By examining the 
equations of motion of the different meson fields (|9.20p . it is clear that the scalar meson field is 
isospin independent. The omega meson field has an implied isospin dependence through the A„- 
coupling to the rho meson field (which is explicitly isospin dependent). Thus the only explicitly 
isospin dependent part of the energy density, defined as e s , is 

1 1 

e s = -^m 2 bl - Av(gvVo) 2 {g p b ) 2 + -g p b (p p - p n ) 

1 f kp 1 [ kn 

H — 77 / dk k 2 V k 2 + to* 2 + I dkk 2 Vk 2 +m* 2 
k 2 Jo vr 2 Jo 

= ~ b o 2 ( m P + 2 A v(gvV ) 2 g 2 p ) + ^g P b (p P - pn) 

l f kp l f kn 

H — 7s / dk k 2 V k 2 + to* 2 + / dk k 2 \J k 2 + m* 2 . 

Jo 7T Jo 

From the equation of motion of the rho meson field (|9.20c|) . bo can be written in terms of t as 

, 1 g P (Pp - Pn) 

2 to 2 + 2A^ 2 (<^o) 2 

1 g P (-t)p 

2 to 2 + 2A^ 2 (^F ) 2 ' 
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The expression for e s /p then becomes 

9 P (-t)p \ 2 (m 2 p + 2A v (g v V ) 2 g 2 \ g p (-t)p (I g P (-t)p 



e 



+ 



p \2 mj + 2A v gj(g v V ) 2 J \ 2p J 2p \2 m 2 + 2A v g 2 (g v V f 

''kn , 1 P fof. 



1 f p 1 f n 

I — 7T- I dk k 2 y k 2 + m* 2 -\ / dk k 2 V k 2 + m* 2 

7T P JO 77 l p Jo 



t 2 p ( 9 2 



8 \m 2 + 2A v g 2 (g v V ) 2 



i r p i r Kn 

+-*- / dkk 2 Vk 2 + m* 2 + dkk 2 Vk 2 + m* 2 . (9.24) 

k 2 P Jo 7T P io 

Thus the defining equation for the symmetry energy coefficient Eq. (|4.ip reduces to, 

1 ( d 2 e s 



° 4 " ~A'dt 2 p) t ^ 



k 3 ( n 2 \ k 2 



+ 



12vr 2 ym 2 + 2A v g 2 (g v V ) 2 J g /^i + m * 2 ' 

after some algebra and using the definition for the baryon density at saturation (|8.25|) . Note 
that the last two terms of Eq. (|9.24p are the same as the two terms in Eq. (|8.26|) and the answer 
can be directly substituted. 

In the Ring convention e s is given by 

\ m l h o + 7,9 P bo{p P - Pn) 



2 p u 2^ 

1 f kp 1 f kn 

-^r / dkk 2 Vk 2 + m* 2 + \ dk k 2 V k 2 + m* 2 , 
Jo J 



and bo by 



1 9p(Pp ~ Pr, 

2 m 2 

1 9p(~t)p 

2 m 2 



The expression for 04 in the Ring convention yields 



a i 



h 3 / \ 2 h 2 
K f 9p\ . K f 



12vr 2 \m p J 6 / k 2 + m * 2 
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9.9 Summary 

In this chapter the equation of state of nuclear matter was calculated using various parameter 
sets of the QHD model of nuclear matter in RMF approximation. 

This information will be used to calculate the equation of state of neutron star matter. 



CHAPTER 10 
Description of neutron star matter 



10.1 Introduction 

Neutron stars are extraordinary objects. Having approximately the same density as atomic 
nuclei but being about 10 20 times larger, these objects have to be approached in a manner 
different from "normal" matter. If the densities of neutron stars are considered, one can con- 
clude that neutron stars could consist, in some part, of normal nuclear matter (protons and 
neutrons). Thus an oversimplified description of a neutron star is that of a giant nucleus with 
mass number of 10 [3]. As mentioned in Chapter Q3 the difference between a neutron star and 
a normal nucleus is that where the nucleus is bound by the nuclear strong force, the neutron 
star is bound by gravity. Thus the description of a neutron star must be anchored in general 
relativity to account for the gravitational effect of not only the mass, but also the energy of 
this ultra-dense object. Using the principles of the special and general theories of relativity the 
Tolman-Oppenheimer-Volkoff (TOV) equation describing a static, spherical symmetric neutron 
star was derived in Chapter [3l To solve the TOV equation, the relationship between the energy 
density and the pressure (equation of state) in the interior of the star must be known. 

The neutron star must also be charge neutral: while the gravitational forces are very strong 
on the large scale, a nett charge would result in very disruptive Coulomb-forces in the neutron 
star [3]. 

The focus in this chapter will be on deriving the equation of state for matter in the neutron 
star interior by considering various physical effects, such as beta-equilibrium as well as the con- 
tribution of the neutron star crust. 

10.2 Crustal effects 

The crust of a neutron star is usually divided into the inner and the outer crust. The outer 
crust ranges from the density where atomic nuclei become fully ionised, at about 10 4 g/cm 3 
(1CP 8 MeV/fm 3 ), up to the density at which the nucleus becomes so neutron-rich that neutrons 
start to drip from the nucleus ( "the neutron drip line" ) [45] . 

The inner crust covers the density range from the neutron drip line up to the density (denoted by 
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Energy density (MeV/fm 3 ) 

Figure 10.1: The BPS equation of state, calculated using the original sequence of 
equilibrium nuclei (labelled as BPS) and the sequence given by the TMA model 
(TMA). 

p c ) at which nuclei cease to exist and nuclear matter undergoes a phase transition to a uniform 
liquid 06]. 

10.2.1 Outer crust 

The BPS equation of state is commonly used to describe the outer crust of neutron stars. 
This equation of state was first described by G. Baym, C. Pethick and P. Sutherland in 1971 
[46] . The departure point of the BPS equation of state is that the nuclei in the crust of a neutron 
star are arranged in a body-centred cubic (bcc) lattice immersed in a sea of free electrons, and 
that the energy of the system can be lowered by the capture of energetic electrons by protons 
(inverse beta decay) [16]. The beta-decay of neutrons in the outer crust is (at this point) blocked 
by the filled electron states [36]. The lattice contributes to lower the energy of the system by 
minimising the Coulomb interaction energy of the nuclei [45]. 

At a specific pressure the nucleus that minimises the chemical potential, /u, of the system is 
deemed the equilibrium nucleus. The pressure range spanned by the outer crust is therefore 
defined by a sequence of equilibrium nuclei. Each equilibrium nucleus defines a certain pressure 
range which is characterised by the maximum mass (energy) density, e max , at which the nucleus 
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can be deemed the equilibrium nucleus. At any given pressure in the outer crust the properties of 
the equilibrium nucleus (either gathered from experimental data or from some theoretical model 
of the nucleus) together with the energy of the electrons determine the energy density. Since the 
pressure in the neutron star crust increases as a smooth function of depth, each transition from 
one equilibrium nucleus to another, will be accompanied by a discontinuity in the energy density 
[46j . This can be clearly seen in Fig. 110.11 where the BPS equations of state represented by the 
sequence of equilibrium nuclei in Tables 110.11 and 110.21 have been plotted. 

The BPS equation of state is described in detail in Refs [H] and [46j . and can be summarised as 
follows: the pressure, P, is given, for an equilibrium nucleus characterised by A (mass number) 
and Z (proton number), by 



• P e is the contribution from the free electrons with fermi momentum (k e ) and mass (m e ) 



P = P e + -wlun 



(10.1) 



where 



• pn is the number density of nuclei, related to the baryon number density (pb) by 



Pb = Ap N 




• wl is the lattice energy per nucleus 



w L = -1.819620 Z 2 e 2 /a, 



where e is the charge of an electron and a is given for a bcc lattice by, 



The energy density (e) is given by 



e = e e + p N (w L + w N ) 



(10.2) 



where 
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• wn is the total energy of an isolated nucleus, and, 

• e e the total electron energy density given by 

e e = dkk 2 ^/k 2 + m e 2 . (10.3) 

Jo 

The equilibrium nucleus is determined by the A, Z configuration that minimizes the chemical 
potential, u } 

_ wn + §u>l + Z^kj + ml 
M= A ' 

The sequence of equilibrium nuclei used in the original BPS paper is shown in Table [T0TT1 where 
the properties of the equilibrium nuclei were calculated in Ref. [35], using the nuclear data set 
that was used in the original BPS paper (Ref. [46] . which was published in 1971, and used nu- 
clear data from the 1960's). Since experimental techniques have undergone some improvement 
since the 1960's, the sequence of equilibrium nuclei between mass (energy) densities of 10 4 g/cm 3 
(1(T 8 MeV/fm 3 ) and 10 10 g/cm 3 (1(T 2 MeV/fm 3 ) can be determined from experimental data 
|45j . The sequence of equilibrium nuclei in Table 110.21 was calculated in Ref. [45] using modern 
experimental nuclear data and nuclear models. The upper part of the table was determined from 
experimental data (compiled by G. Audi et ah), while the lower part was calculated using a 
modern relativistic nuclear field theory model by L.S. Geng et al. called TMA (see Ref. [45] and 
references therein). 

The BPS equation of state can be determined from either Table 110.11 or Table 110.21 It can 
be done by either interpolating the energy density at a given pressure, or at a specific pressure 
(and hence equilibrium nucleus) the contribution of the electrons to the pressure can be deter- 
mined from Eq. (jlO.ip and the fermi momenta of the electrons (k e ) obtained. If k e is known, e e 
can be determined from Eq. ()10.3p . Using the value of e e and the properties of the equilibrium 
nuclei, the energy density can then be obtained from Eq. (|10.2p . 

10.2.2 Inner crust 

The BBP equation of state (by G. Baym, H. A. Bethe and C. J. Pethick [37]) assumes that 
at densities above the neutron drip line, up to p c , nuclear matter consists of a lattice of nuclei 
immersed in a gas of free neutrons. The equation of state is calculated by matching the pressure 
and the energy density of the lattice nuclei to that of a free neutron gas [47J. The equation of 
state is given in Table [TITH (at the end of the chapter). 
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30 
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50 
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26 
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82 
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25.25 
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2.23 x 10 35 


120 Sr 


38 


82 


939.57 


26.19 


4.32 x 10 11 


7.87 x 10 29 


2.59 x 10 35 


118 Kr 


36 


82 



Table 10.1: Sequence of equilibrium nuclei, from Ref. |45j, that are used to reproduce 
the original equation of state of the outer crust of a neutron star as in Ref. |46j. 
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Table 10.2: Sequence of equilibrium nuclei, calculated by using modern atomic data 
and the relativistic TMA model of the nucleus, from Ref. |45j, that are used to 
calculate the equation of state of the outer crust of a neutron star. 



The value of the density at which the inner crust becomes a uniform liquid (p c ) depends on 
the model (and parameter set) used to describe the matter of the neutron star crust and interior. 
J. Carriere et al. calculated the value of p c for different QHD parameter sets in Ref. |48| . In 
Ref. [H] the authors used the relativistic random-phase approximation to analyse if and when a 
system consisting of protons, neutrons and electron would undergo a second-order phase transi- 
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tion. For the inner crust they assumed a polytropic equation of state, 



P(e) = A + Be*/ 3 



(10.4) 



where A and B are constants that are determined so that the above equation of state matches 
the equation of state for the outer crust (pressure and energy density) at the neutron drip line 
as well as the equation of state of the liquid interior at p c . 

10.3 Equilibrium conditions 

Free neutrons are not stable particles and have a half- life of about 10 minutes [59]. The 
neutron will decay into a proton, electron and an antineutrino (beta-decay): 



Therefore free neutrons in the neutron star will undergo beta-decay creating a proton and neu- 
tron. Since a neutron stars exists for a period longer than the half-life of a neutron, a significant 
number of neutrons in the star will undergo beta-decay. Since a neutron star consists of a huge 
amount of neutrons, the number of electrons (resulting from the beta-decay of neutrons) present 
in the star will also be quite significant. Due to the densities of the electrons, some electrons 
will attain relativistic energies, making the inverse process (capture of a electron by a proton to 
form a neutron) energetically favourable, i.e. 



The two decay processes mentioned above represent two competing processes and since the 
neutron star has a lifetime that is long compared to the timescale of these decay processes, the 
matter in the star has time to reach beta-equilibrium. The beta-equilibrated state is the state 
where, for a fixed baryon number density (p = p p + p n ), the proton and neutron number densities 
(p p and p n ) are such that the energy density (e) of the system is at a minimum. From Ref. [4J, 
the condition for beta-equilibrium can be expressed in terms of the chemical potential (Fermi 
energy), p. If p x refers to the chemical potential of a specific particle specie, x, then 



The effects of the neutrinos are not considered in the beta-equilibrated state, for their mean-free 
path is longer than the radius of the star [3J, and therefore the beta-equilibrated state is described 



n — > p + e + v e . 



p + e — > n + v e . 



de 



'X r, 

OPx 
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by 

= fJ-p + Me • (10.5) 

The number density of particle x (p x ) can be expressed in terms of the fermi momentum (k x ) as 

7,3 

_ K x 

Px - 3vr2 ' 

and from the above k x is 

k x =7T 2 (37^)" 2/3 , 

which is a handy result when evaluating the chemical potentials of different particles. 

Using the expression for the energy density (|10.11a|) . the chemical potentials of the different 
particles can be expressed as 

[i n = \Jk 2 + m* 2 + g v V - -g p b 



1 



= \J k l + m * 2 + 9vVq + -g p b 



The condition for beta-equilibrium (|10.5|) can therefore be generalised, for all parameter sets 
considered in this work, to 



k 2 + m* 2 = J k 2 + m* 2 + g p b + yjk 2 + m e 2 . (10.7) 
In the Ring convention, the expression for beta-equilibrium is 



k 2 + m* 2 = Jkj+ m* 2 - 2g p b + \Jk 2 e + m e 2 . (10. 



In ultra-dense matter electrons can attain ultra-relativistic energies and therefore it may become 
energetically more favourable to populate muon states. Muons have the same charge as an 
electron, but a mass of 105.7 MeV |50| . In this work it is assumed that muon states will be 
populated once the electrons reach energies comparable to the mass of the muon. Muons are 
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unstable particles and undergoes the following decay reaction [4]: 



e +v e + v ll . 



Equilibrium with regards to the above reaction implies that 



Me = /V 



(10.9) 



where the muon chemical potential is given by 



Since the neutron star is assumed to be neutral, the number of protons must be equal to the 
number of electrons present. In terms of the fermi momenta of the species this means that 



but if muon-states are also populated, the condition for charge equilibrium will have to be 
expressed as 



Neutron stars are in general equilibrium when the conditions of (|10.5p . f|10.9[) and (jlO.lOp are 
satisfied. 

10.4 Equation of state 
10.4.1 Without a crust 

The energy density and the pressure of neutron star matter in general equilibrium, assuming 
a QHD description for nuclear matter that contains protons, neutrons, electrons, muons and 



Pfi + Pe = Pp 



which implies that 




(10.10) 
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Figure 10.2: Behaviour of the equation of state (EoS) of the neutron star interior 
when different equations of state for the inner crust are used. 

mesons (a, ui and p) is, from (|9.22[) : 



+ 



1 A 1 C 

2 m2 s4 + |(s^o) 3 + ^(^o) 4 - - 2 mlV% - ±(g v V ) 4 



7, m l b o ~ A v{gvV ) 2 (g p b ) 2 + g v V (p n + p P ) + ^g P bo{p P ~ Pr, 
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7T ./ 



dk k 2 \/k 2 + m* 2 



+ s^y dkk 2 ^[w+ 



(10.11a) 



and, 



2 ±2 

m s (p 



^(9sM 3 ~ ^( 9 sM 4 + \mlV 2 + ^{g v V f 



+ ^m 2 p b 2 + A v (g v V ) 2 (g p b ) 2 



+ 



1 

3^ 



k 



1 



+ £ 



dk — H ^ / dfc 

Vk 2 + m* 2 3^ J 

1 f^A 



37r 2 y 



(10.11b) 



where A denotes the lepton (electrons and muons) state. 
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10.4.2 Including a crust 

In this work the BPS equation of state is used for the outer crust. Both the original sequence 
of nuclei (of Table [TOTi]) and modern calculations of the sequence of nuclei using the TMA model 
(as in Table 110, 2p are used. The difference between the resulting equations of state is shown in 
Fig. 110.11 When the values of Table 110.11 are used for the sequence of equilibrium nuclei, it will 
be denoted by BPS. When those of Table [10721 are used, it will be denoted by TMA. 



For the equation of state of the inner crust the BBP equation of state was initially used, but 
since only the energy density or the pressure can be matched at p c the equation of state has a 
rather large discontinuity, as can be seen in Fig. 110.21 As such the BBP equation of state was 
not used in further calculations presented in this work. Rather the polytropic equation of state 
(jlO.ip for the inner crust is used, since the values of A and B can be calculated to the equation 
of state of the liquid interior at p c . Due to time constraints, the calculation of p c for the different 
QHD parameter sets was not repeated. Rather the published values for p c , where available, were 
used. For the parameter sets for which values of p c are not available, the value of the pressure 
at p c was estimated to be 4.982 x 10 32 dyne/cm 2 (0.311 MeV/fm 3 ). This pressure is simply the 
average pressure, corresponding to the published values for p c , in the specific parameter sets. 
The published values for p c , together with the corresponding pressures, are listed in Table [T0~3l 



When crustal effects are included, the equation of state of the liquid interior of the star is 
given by the energy density and the pressure derived for the QHD description of the neutron 
star interior (jlO.lip . 

The equation of state of the neutron star interior, when crustal effects are included, can therefore 
be summarised as 



P(e) 



^BPs(e) 



if e min — e — e outer 



A + Be 4 / 3 if e uter < e < e c 
-Pqhd(c) if e > e c 



10.12) 



10.5 Neutron star constraints 

As with any investigation into natural phenomena, experimental and observational inputs 
are crucial to the development of the description of a physical system. In the case of neutron 
stars these inputs could stem from experiments or astrophysical observations. Recently a testing 
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Parameter set 


Pc 


Pressure at p c 




(in MeV-fm- 3 ) 


(in dyne-cm -2 ) 


FSUGold 


0.076 


6.56 x 10 s2 


NL3 


0.052 


3.42 x 10 32 



Table 10.3: The values for transition (number) density from the inner crust to 
the liquid interior of the FSUGold and NL3 parameter sets, from Ref. |44|. The 
corresponding pressures (at p c ) were calculated in this work, assuming that the 
neutron star is in general equilibrium. 

scheme based on various different constraints derived from heavy-ion collisions and astrophysical 
observations has been proposed by T. Klahn et al. [51] • In this section the constraints proposed 
in Ref. [51] and other publications, applicable to the scope of this work, will be discussed. 

10.5.1 General theory of relativity 

From the Tolman-Oppenheimer-Volkoff equation (|3.1|) it is implied that the radius of the star 
(R) cannot be smaller or equal to 2GM. If this condition is not obeyed, the change in the 
pressure is infinite or increases as one moves from the centre of the star to the boundary. 2GM 
defines a massive object's Schwarzschild radius and if the radius of the object lies within its 
Schwarzschild radius, i.e. 

R < 2GM, 

then the object has collapsed under its gravity and is a black hole [I]. 

10.5.2 Observational constraints 

The mass of a neutron star in a binary system can be quite accurately measured [5] and the 
mass of various neutron stars have been determined. The list of known masses of neutron stars 
is given in Fig. 110.31 

The known neutron star masses cover a narrow range, between about 1.4 and 1.5 M & (solar 
mass). This narrow mass range may be due to special circumstances that govern the evolution 
of binaries, but no thorough explanation has been put forward to date [5j. The heaviest known 
neutron star is PSR J0751+1807 with a mass of 2.1 ± O.2M |52j . The reproduction of this mass 
by any model for the neutron star interior (equation of state) has been hailed as a very important 
constraint in determining whether certain models for the neutron star interior are applicable or 
not |51| . J. M. Lattimer and M. Prakash have however warned in Ref. [5] that caution should be 
exercised before the currently known heavy neutron stars are taken as firm evidence that such 
heavy neutron stars exists, due to uncertainties in the measurements. 
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Figure 10.3: Table of masses (with error bars) of known neutron stars. Taken from 
Ref. [5]. 



To determine the radius of a neutron star is much more difficult than the mass and there are large 
uncertainties in the measurements [3|. To date, the only neutron star whose mass and radius has 
been determined is EXO 0748-676 [53]. A. R. Villarreal and T. E. Strohmayer established that 
the mass of EXO 0748-676 must lie between 1.5 and 2.3 M & while the radius could lie between 
9.5 and 15 km [53]. F. Ozel determined the lower limits of the mass and radius of EXO 0748-676 
to be 2.10 ± 0.28 M and 13.8 ±1.8 km [53]. 



These considerations are used in this work to evaluate different equations of state for the neutron 
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star interior. 

10.6 Summary 

In this chapter a description of a neutron star was presented, based on certain assumptions 
about the crust of the neutron star as well as equilibrium conditions in the interior of the star. 
The equation of state of neutron star matter was also derived. 
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Table 10.4: The BBP equation of state of matter in the inner crust of a neutron 
star, as calculated in Ref. [47] . 



CHAPTER 11 
Results and discussion 



The results in this chapter were generated by writing a computer program to solve the applicable 
equations in the FORTRAN90 programming language. A description of the programs as well as 
the actual equations that were used in the computations, are given in Appendix [Al 

11.1 Properties of saturated nuclear matter 

The properties of saturated nuclear matter were calculated using the different parameter sets 
discussed in Chapters [8] and KQHD-I, NL-SH, TM1, TM2, NL3, PK1 and FSUGold). The cal- 
culated values are compared to the published ones for the different parameter sets. The results 
are listed in Tables 111.21 and 111.31 The properties of saturated nuclear matter calculated in 
this work shows good agreement with the published values. 



Parameter set 




Psat(fm a ) 






E fe (MeV) 






Published 


Calculated 


Observed 


Published 


Calculated 


Observed 








0.153 






-16.3 


QHD1 [22] 


0.148 


0.148 




-15.75 


-15.75 




NL-SH @D] 


0.146 


0.146 




-16.328 


-16.346 




TM1 [H] 


0.145 


0.145 




-16.3 


-16.3 




TM2 [H] 


0.132 


0.132 




-16.2 


-16.2 




NL3 02] 


0.148 


0.148 




-16.299 


-16.240 




PK1 [43] 


0.148195 


0.148192 




-16.268 


-16.268 




NL3 [55] 


0.148 


0.148 




-16.24 


-16.24 




FSUGold [551 


0.148 


0.148 




-16.30 


-16.28 





Table 11.1: Comparison between the published values and the ones calculated in 
this work for the saturation density of nuclear matter and the binding energy per 
nucleon at the saturation density. The references after the name of the parameter 
set indicate the source of the published values. 



To determine the accuracy with which the FORTRAN90 program calculates the properties of 
saturated nuclear matter, these properties were also calculated using the NL-SH, TM1 and TM2 
parameter sets. However, from this point onwards only the QHD-I, NL3, PK1 and FSUGold 
parameter sets will be considered. 

The binding energies around the saturation density (at a fermi momentum of 1.30 fm _1 ) of 
the QHD-I, NL3, PK1 and FSUGold parameter sets are shown in Fig. 111.11 From this plot it 
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Parameter set 


Published 


K(MeV) 
Calculated 


Observed 


Published 


a 4 (MeV) 
Calculated 


Observed 








234 






32.5 


QHD1 [22] 


not available 


546.55 




uot available 


19.3 




NL-SH 0D] 


354.95 


355.34 




36.1 


36.1 




TM1 [31] 


281 


281 




36.9 


36.9 




TM2 [31] 


344 


344 




35.8 


36.0 




NL3 [32] 


271.76 


271.52 




37.4 


37.4 




PK1 [43] 


282.644 


282.685 




37.641 


37.640 




NL3 [55] 


271 


272 




37.3 


37.3 




FSUGold [55J 


230 


230 




32.6 


32.6 





Table 11.2: Comparison between the published values and the ones calculated in 
this work for the compressibility and symmetry energy of nuclear matter at satura- 
tion. The references after the name of the parameter set indicate the source of the 
published values. 



Parameter set 


m*/M n 


m*/M n 




Published 


Calculated 


Published 


Calculated 


QHD1 [22] 


0.54 


0.54 


0.54 


0.54 


NL-SH [30] 


0.60 


0.60 


0.60 


0.60 


TM1 [31] 


0.634 


0.634 


0.634 


0.634 


TM2 [31] 


0.571 


0.571 


0.571 


0.571 


NL3 [32] 


0.60 


0.60 


0.60 


0.60 


PK1 [43] 


0.605525 


0.605526 


0.604981 


0.604983 


NL3 


not available 


0.60 


not available 


0.60 


FSUGold 


not available 


0.61 


not available 


0.61 



Table 11.3: A comparison of the published values and the ones calculated in this 
work of the ration of the reduced (baryon) mass to the baryon mass at satura- 
tion. The references after the name of the parameter set indicate the source of the 
published values. 

is clear that the different parameters sets have very different high density behaviour. All the 
parameter sets give almost exactly the same values for the binding energies below saturation, 
but above saturation the values diverge. (The QHD1 parameter set was fitted to reproduce a 
binding energy at saturation of -15.75 MeV and therefore its description differs slightly from that 
given by the other parameter sets [22]-) These discrepancies in the description of the binding 
energy at high densities emphasise that more information about dense matter is needed to better 
constrain the descriptions of nuclear matter. 

11.2 Neutron star properties 

In this work the NL3, PK1 and FSUGold parameter sets will be used to compare different 
calculated neutron star properties. The FSUGold parameter set will be the default. 
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k p (fm- 1 ) 

Figure 11.1: A comparison of the binding energy per nucleon of nuclear matter of 
different parameter sets. 

11.2.1 Inclusion of different particles 

As a first approximation the neutron star can be considered to be exclusively composed out 
of neutrons. This satisfies the condition of charge neutrality, but not the condition of beta- 
equilibrium, as expounded upon in Chapter [TUJ By including the effects of beta-equilibrium the 
maximum mass of the neutron star sequence is reduced, i.e. the equation of state is softened. As 
was explained in Sec. l4.1.-H a hard equation of state is when the energy density increases rapidly 
with an increase in pressure and a soft equation of state is when the increase in the energy density 
is more gradual. The equation of state is softened by the inclusion of different particles in the 
description of neutron star matter, because if it becomes energetically favourable to populate the 
states of, for example muons, the lowest energy muon states will be populated first, instead of 
the electron states, thereby reducing the increase in the energy density as the pressure increases. 
Fig. 1 1 1 . 2 1 is an illustration of the softening of the equation of state by the inclusion of additional 
particles to the description of neutron star matter. If the neutron star is further considered to 
not only be in beta-equilibrium, but in general equilibrium as well (thus also including muons, 
instead of just protons, neutrons and electrons), the maximum mass of a particular sequence is 
further reduced. 
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Figure 11.2: Softening of the equation of state due to the inclusion of different 
particles to the description of the neutron star interior. 



Parameter set 


Maximum mass (M©) 




Pure neutron matter Beta-equilibrated matter Matter in general equilibrium 
(only neutrons) (+ protons and electrons) (+ muons) 


NL3 
PKl 

FSUGold 


2.902 2.786 2.769 
2.546 2.342 2.308 
1.800 1.737 1.718 



Table 11.4: Reduction of the maximum neutron star mass as different particles are 
included in the description of the neutron star interior. 



Table 111.41 lists the maximum masses of a pure neutron neutron star, one in beta-equilibrium 
and a neutron star in general equilibrium, for different parameter sets. Fig. 111.31 shows the 
effect of the softening of the equation of state on the mass-radius relationship (equation of state 
described by the FSUGold parameter set). 

11.2.2 Crustal effects 

The inclusion of crustal effects does not affect the mass of the neutron star, when compared 
to the case where crustal effects have not been included. However, the radii of neutron stars are 
in general larger when crustal effects are considered. 

The reason that the inclusion of crustal effects do not influence the mass of the neutron star 
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Figure 11.3: Influence, due to the inclusion of various particles in the description 
of the neutron star interior, on the mass-radius relationship of neutron stars. 



Parameter set 


TMA 


BPS 




Maximum mass (Mq) 


Radius (km) 


Maximum mass (M©) 


Radius (km) 


NL3 


2.768 


13.337 


2.768 


13.331 


PK1 


2.308 


12.761 


2.308 


12.752 


FSUGold 


1.718 


10.893 


1.718 


10.863 



Table 11.5: Maximum masses and radii of different neutron star sequences. (The 
neutron star matter is considered to be in general equilibrium and crustal effects 
are included in the description.) 



is that it does not affect the central pressure in the neutron star. The central pressure corre- 
sponds to the amount of overlying mass. The effect of the inclusion of crustal effects is rather 
to lower the pressure gradient near the boundary of the star, thereby enlarging the radius of the 
star. Fig. 111.41 shows the neutron star mass-radius sequence, for an equation of state based on 
the FSUGold parameter set, where different descriptions for the outer crust (TMA and BPS), as 
well as the polytropic equation of state for the inner crust have been included. Fig. 111.51 shows 
the plot of the neutron star mass against the central density of this calculation. 
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Figure 11.4: Sequence of neutron star mass-radius relationships, where outer and 
inner crustal effects have been included in the description of the neutron star in- 
terior. The TMA plot refers to the BPS equation of state where the sequence of 
equilibrium nuclei is given by Table [T0.2L while the BPS plot relates to that of Table 

[TEE 
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Figure 11.5: Plot of neutron star mass versus central density, where crustal effects 
have been included, corresponding to Fig. 111.41 
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Figure 11.6: Sequences of neutron star mass-radius relationships corresponding to 
equations of state based on different parameter sets. All equations of state shown 
include a description of the inner crust (polytropic equation of state, matched at the 
appropriate densities for the transition to the liquid interior) and the BPS equation 
of state (specifically the TMA rendition) for the outer crust. The liquid interior of 
the star is assumed to consist of protons, neutrons, electrons and muons in general 
equilibrium. The horizontal (pink) lines correspond to the constraints imposed by 
the mass measurement of PSR J0751+1807. The diagonal dotted (purple) line 
represents the possible mass-radius relationships of EXO 0748-676, as calculated by 
F. Ozel |53j, while the diagonal solid (purple) line corresponds to the calculation 
made by A. R. Villarreal and T. E. Strohmayer |54j . The diagonal solid (black) line 
in the upper left-hand corner denotes the constraint imposed on the mass-radius 
relationship of a neutron star by the general theory of relativity. 

II. 2.3 Observational constraints 

As argued in Section 110,5.21 observational evidence and theoretical constraints have to be 
taken into account to determine the appropriate descriptions of the neutron star interior. Fig. 

III. 61 shows the plot of mass-radius sequences corresponding to different equations of state. 

From Fig. 111.61 it can be seen that the NL3 and PK1 parameter sets provide a description 
which is in good agreement with current observational data, while the FSUGold description may 
seem to be too soft, if the maximum masses of each sequence (as given in Table I11.5P are con- 
sidered. 
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J. Piekarewicz has pointed out these apparent deficiencies of the FSUGold parameter set in 
Ref. [55], but also stated that it would be rather naive to only consider neutron star constraints 
to determine the equations of state of dense matter. The difference between NL3 and PK1 
parameter sets, is that the PK1 set includes an additional self-coupling in the omega meson 
field (V^). Comparing the PK1 and FSUGold parameter sets, it will be seen that the FSUGold 
parameter set contains a further coupling between the two vector boson fields (b^ and V^), A„, 
which was introduced to soften the symmetry at high density [44]. However, the density depen- 
dence of the symmetry energy is currently unknown, so A„ could not yet be firmly constrained 
[55] . However, the Parity Radius Experiment (PREX) aims to accurately measure the neutron 
radius of 208 Pb, which would provide an accurate determination of the density dependence of the 
symmetry energy [3]. PREX will be conducted at Jefferson Laboratories in the United States 
and is currently scheduled for February 2009 [56J. 
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Conclusions 



In this work the equations of state of neutron star matter were derived using the quantum 
hadrodynamics model of nuclei and nuclear matter by the application of the relativisitc mean- 
field approximation for the NL3, PK1 and FSUGold parameter sets. It was also shown (as was 
done in similar works, such as Refs |48[ [55] and [57]) that these derived equations of state can 
be applied to calculate the mass-radius relationship of a neutron star and that these calculated 
values show good agreement with current observational data. However, further measurements 
of the properties of nuclear matter as well as that of neutron stars are necessary to obtain a 
complete description of dense nuclear matter. 

As argued by various authors, a description of dense matter should not be based on a single 
set of observables or a class of observables, but rather the most applicable description of dense 
matter would be the one that describes a plethora of different observables accurately |51[ [55] . Un- 
fortunately our knowledge of dense matter, in whatever form, is currently limited. Fortunately it 
is foreseen that with upcoming experiments at current and new particle accelerator facilities, as 
well as the development of new telescopes, particularly radio-telescopes, much will be uncovered 
and discovered about terrestrial and extra-terrestrial dense matter in the near future. To gather 
a more complete picture of dense matter, works, such as this one, would be increasingly im- 
portant in providing a link between laboratory physics and astrophysics. However, through the 
study of nuclear physics alone, this picture would not be easily attained and therefore the need 
for cooperation between astro- and nuclear physicists is greater than ever. 

The application of relativistic mean-field theory to the quantum hadrodynamics has therefore 
shown itself to be relevant to the description of finite nuclei and nuclear matter [22], as well as 
the description of neutron stars, based solely on a handful of coupling constants. 
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APPENDIX A 
Code Documentation 



All numerical problems in this work were solved using programs that were written for this 
purpose in FORTRAN90 programming language. This appendix will give a brief overview of 
these computer programs. 

A.l Introduction 

To obtain the mass-radius relationship of a neutron star using a specific parameter set in QHD 
is quite a complex physical and mathematical problem. Even after using RMF to approximate 
the description of nuclear matter, the equations that have to be solved are rarely exact, at most 
self-consisted. The main problem of applying a description of nuclear matter to neutron stars 
can be split into the following components: 

1. obtaining the values of the unknown meson fields in the RMF approximation, 

2. calculating the equation of state of neutron star matter, and, 

3. using the equation of state to solve the TOV equation. 

In the following sections different parts of the code pertaining to the above-mentioned problem 
will be discussed. Whenever a function or subroutine is mentioned, it will be written in bold 
and subroutines are denoted by the prefix sub- and functions by the prefix func_. 

A. 2 sub_fields 

Sub_fields is the subroutine used to calculate the RMF values of the meson fields (cf>o, Vq 
and bo). 

A schematic representation of this subroutine is given in Fig. IA.1I As imput sub_fields receives 
values of the proton and neutron fermi momenta, k p (denoted by k_p) and k n (k_n), and the 
choice of the parameter set (para_choice) for which the RMF meson field values have to calcu- 
lated. For each specific value of k_p and k_n sub_conversion is called. Sub_conversion assigns 
all the coupling constants and masses as well as the baryon number density (density _B) as global 
parameters for the specific choice of para_choice. Once this is done the value of the RMF meson 
fields can be calculated. This is done by solving the self-consistent expressions for the different 
meson fields (IQ0|) . 
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Receives 

proton fermi momentum (k_p) 
neutron fermi momentum (k_n) 
parameter set choice (p_c) 



Returns 
x_p & x_n 



Assigns 

RMF meson fields 
sigma_0* 
omega_0* 
rho 0* 



k_p 
k_n 



sub conversion 



Ring 



Walecka 



k_p 
k_n 
p_c 



coupling constants* 
and 

masses* 
algo_choice* 



Solve self-consistent eqns. 
for RMF meson fields by 
incremental search 



x_p 
x_n 

sigma_0* 
omega_0* 
rho 0* 



Figure A.l: Illustration of the subroutine sub_fields (global variables are denoted 
by an asterik). 



The values of the Vq and &o are determined by comparing the left- and right-hand sides of 
Eqs (|9.20bj) and (|9.20cp . The value of the field on the left-hand side of the equation is incremen- 
tally increased and then substituted into the expression on the right and the values on the left 
and right are compared. The RMF meson fields get assigned once the difference between the left 
and right sides are acceptability small (about 8.0 x 10~ 7 fm -4 ). 

In the case of (pQ the value of the reduced mass (m*) is calculated, from which </>o is then 
calculated using Eq. ()9.7p . An expression for m* in terms of <f)Q can be easily deduced from the 
expressions for 4>q (|9.20ap and m* (|9. T[) and yields 




(A.l) 



The integral in the above expression is evaluated using the analytical expression for the integral 
given in Appendix [Bl 

After checking one final time if the self-consistent equations are satisfied to the desired accuracy 
(this is necessary, especially since in the FSUGold parameter set there is a coupling between Vq 
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Receives choices for 
ciustal effects (crust) 
parameter set (para) 
equilibrium conditions (eqm) 



crust 

para 
eqm 




Figure A. 2: Illustration of the subroutine sub_epNS (global variables are denoted 
by an asterik). 



and 60 which complicates the calculation) the values of 4>q [denoted by sigma_0 (sic) throughout 
the code], m* (m_reduced), Vq (omega_0) and 6q (rho_0) are assigned (as global variables). 



kp 
m* 



and 



kp 
m* 



are also assigned. x p and x n are used in the calculation of the energy density and the pressure 
and stem from the analytical solution to the integrals Appendix [B] in the expression for the 
energy density and the pressure (|9.22|) . 



The differences between the Walecka and Ring conventions (see Sec. I9.3P necessitate the defini- 
tion of a different subroutine to calculate the RMF meson field values in the Ring convention. 
Sub_fieldsl is very similar to sub_fields, but one of the few differences is that a numerical 
approximation is made to calculate the integral in the expression for m*. 

A.3 sub_epNS 

This subroutine calculates the arrays containing corresponding values of the energy density 
and the pressure to be used in sub_TOV, which solves the TOV equation. Sub_epNS is 
represented graphically in Fig. IA.21 
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Specific choices of the crust (crust_choice), equilibrium (beta of general, eqm_choice) and the 
parameter set (para_choice) serve as input to sub_epNS. 

For a specific choice of the crust sub_pcutoff determines the transition between the inner crust 
and the liquid interior and/or between the outer and inner crust. It returns a value of k n and k p 
(k_cutn and k_cutp) at the corresponding pressure or number density of the transition. 

If the description of matter in general equilibrium is required, sub_muon returns the value 
of k n (k_max_no_muons) above which the muons states will be populated. Sub_muon also as- 
signs the minimum values of k n (k_n_beta_min) at which the matter can be in beta-equilibrium. 
The value of k_n_beta_min is not calculated by sub_muon, but was determined by hand and 
then hard coded to be 1.49 x 10 _1 fm _1 . Sub_beta is used to determine the values of the proton, 
electron and muon fermi momenta (kp, ke and kmu) corresponding to a neutron fermi momen- 
tum, kn (which of course this calculation is parameter set specific). If only beta-equilibrated 
matter is desired, sub_beta calls sub_beta2, which does the necessary. 

If no crustal effects are considered, sub_fields is called to determine the values of 4>q (sigma_0), 
loq (omega_0), 60 (rho_0), x p (x_p) and x n (x_n). These values, along with the values of k e (k_e) 
and k^ (k_mu), are then used by func_epsNS and func_presNS to calculate the energy and 
the pressure of neutron star matter, using Eqs (I9.22a|) and (I9.22b[) respectively. 

The equation of state for the outer crust is calculated by func_bps. For the choice of the 
sequence of equilibrium nuclei (either TMA from Table [HOI or BPS from Table [IOTT]) and specific 
value of the pressure in the outer crust, func_bps returns the energy density. For the inner crust 
the BBP equation of state can be used or the polytropic equation of state [denoted by RPA 
in the program and given by Eq. (|10.4p ]. The BBP equation of state is calculated by reading 
the values given in Table 110.41 from file (BBPE.txt and BBPP.txt for the energy density and 
the pressure respectively) and determined by interpolating the energy density at a specific pres- 
sure. The polytropic equation of state is calculated by first determining the values for A and B. 
Using equation 110.41 the energy density is calculated for each input value of the pressure. The 
equation of state of the liquid interior of the neutron star is given by calling func_epsNS and 
func_presNS for values of k n (kn) greater than kcutn. 
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The output of sub_epNS is a file, ep.txt, which contains the pressure in the first column and 
the energy density in the second. The following data is also printed to file: The value of k n (kn) 
and k p (kp) to kn_kp_ep.txt; k e (ke), k^ (kmu) and density _B (p p + p n ) to ke_kmu_ep.txt, while 
fields_ep.txt contains the values of <f>o (sigma_0), loq (omega_0) and 60 (rho_0). 

A.4 sub_TOV 

This subroutine solves the TOV equation. As input sub_TOV receives the neutron fermi 
momentum (k_fn) assumed to describe the matter at the centre of the star, as well as the choices 
for the parameter set, crust and equilibrium conditions. The equation of state have to be read 
into the pres and eps arrays as global variables, this has to be done before sub_TOV is called. 



Sub_TOV solves the TOV equation using the Runge-Kutta method to solve a coupled dif- 
ferential equation. 

For a given coupled differential equation, such as Eqs (|3.ip and ([3. 2ft . 



and 



dP 

dr 



dM 
dr 



F(r, P, M) 



GeM 



1 + 



G(r, P, M) 



1 + 



4irr 3 P 
M 



2GM 



(A.2a) 



(A.2b) 



the Runge-Kutta method defines a single integration step of size h by [59 



P n+1 = P n + \k\ + \k 2 + \k% + -U4 
boob 

M n+l = M n + -d 1 + -d 2 + -d z + -d i 
boob 
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with 

ki = hF(r n ,P n ,M n ) 

k 2 = /iF(r n + ty2,Pn + fciAM n + di/2) 

k 3 = hF(r n + h/2,P n + k 2 /2,M n + d 2 /2) 
k 4 = hF(r n + h,P n + k 3 ,M n + d 3 ) 

di = hG(r n ,P n ,M n ) 

d 2 = hG(r n + h/2,P n + k 1 /2,M n + d 1 /2) 

d 3 = hG(r n + h/2,P n + k 2 /2,M n + d 2 /2) 

di = hG(r n + h,P n + k 3 ,M n + d 3 ). 

If crustal effects are taking into consideration sub_TOV uses the above method to integrate 
the TOV equation until the pressure [-P(r)] becomes 10 -8 MeV/fm 3 , which denotes the value 
of the pressure were the crust of a neutron star starts, from Ref. [15]. If crustal effects are 
not considered sub_TOV integrates until the pressure [P(r)] becomes 2.0 x 1(T 5 MeV/fm 3 (the 
lowest pressure at which matter can be in beta-equilibrium). 

Due to the units of the variables care should be taken with the expressions for F(r, P, M) and 
G(r, P, M) as pointed out in Refs |20|. 158] . M and r are of the order of a couple solar masses 
and kilometres respectively, while P and e are in MeV/fm 3 . Thus to be able to handle these 
equations numerically, they have to be scaled. An appropriate choice would be to describe M in 
terms of solar masses, i.e. 

M = MqM 

where 

• Mq is a solar mass (in MeV) as given in Table [TT| and, 

• M is a dimensionless quantity 

as well as to express r in units of metres. 

Defining the following quantities 

• G = G' ■ 10 -30 with G the gravitational constant in MeV/fm 3 from Table [TTT| and, 
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• M & = M' & - 10 45 , 
(1A.2|) can be rewritten as 

dP 

dr 



and 



dM' 
dr 



F(r, P, M) 



G'eM'M' & 



G(r, P, M) 
Aver 2 



p 

e 



1 + 



4nr 6 P 
M'qM' 



2G'M' Q M 



-i 



(A.3a) 



(A.3b) 



The above equations are the ones used by sub_TOV to calculate the mass-radius relationships 
of neutron stars in this work. 



A. 5 List of modules 

For completeness sake the modules in the FORTRAN90 program will be listed. (If the sub- 
routine or function has not been mentioned previously a quick explanation will be given.) 



The module mod_comp_val_RING contains 

• sub_fieldsl, 

• func_rhosI , which calculates the integral in the expression for 0o numerically, and, 

• func_intrhos calculates the integrand of the numerical integration in func_rhosI, while, 

• func_sigmal, 

• func_omegaI, and, 

• func_rhoI are the self-consistent expressions for the meson fields called by sub_fields. 
The module mod_comp_val_WALECKA contains 

• sub_fields, 

• func_rho_s, which calculates the integral in the expression for 0o numerically, and, 

• func_int_rho the integrand of the numerical integration in func_rho_s, while, 

• func_sigma, 
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• func_omega, 

• func_rho are the self-consistent expressions for the meson fields called by sub_fields, and, 

• print .parameters that prints the value of the coupling constants and masses of a particular 
parameter set to the screen. 

The two modules mentioned above contains all the subroutines and functions used in calculating 
the RMF values of the meson fields and refers to the two different conventions used in this work. 

The module mod_EoSFSU contains 

• func_epsNS, 

• func_presNS, 

• func_epsFSU, 

• func_presFSU, 

• func_inteps which contains the integrand for the numerical integration of the integral in 
the expression for the energy density (this function was not used), 

• funcJbps, and, 

• func_getke which calculates the electron fermi momentum for a given pressure in func_bps; 
whereas the module mod_gauleg contains 

• sub_gauleg, which calculates the roots and weights for the numerical integration, 
and mod_splint contains the interpolation functions 

• sub_csplin which calculates the coefficients for the numerical interpolation, and, 

• func_cseval which does the actual interpolation using the coefficients generated by sub_csplin. 
The module mod_globals contains all the global parameters as well as 

• sub_choice a simple subroutine that prints the choices made to the screen, and, 

• sub_conversion, 
while mod_obs contains 
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• sub_sat, which calculates the properties of saturated nuclear matter for a choice of parameter 
set and prints the results to the screen and to the file sat_prop.txt, 

• sub_symm_ep calculates the equation of state of symmetric nuclear matter, used in 
sub_sat, and, 

• sub_epNS. 

mod_procedures houses all kinds of functions and subroutines used by other modules, such as 

• subJbeta, 

• sub_beta2, 

• sub_muon, and, 

• sub_pcutoff, 

which were all introduced previously. 

Finally mod_TOVFSU contains 

• sub_TOV, 

• func_TOVlFus, and, 

• func_TOVlGus which contains the expressions for the right hand side of the equations 
used to solve the TOV equation (|A,3j) . 

A. 6 Codes for different choices 

The codes for choices of the parameters sets, crustal effects and equilibrium condition are: 
Parameter sets 

• 1 FSUGold 

• 2 QHD1 

• 3 NL3 (Walecka convention) 

• 31 NL3 (Ring convention) 

• 51 PK1 

• 61 NL-SH 
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• 71 TM1 

• 81 TM2 
Crustal choices 

• if no crustal effects included 

• 1 if only outer crust, BPS equation of state based on TMA sequence of equilibrium nuclei, 

• 2 if only outer crust, BPS equation of state based on BPS sequence of equilibrium nuclei, 

• 31 if the BPS (TMA) equation of state for outer crust and BBP equation of state for the 
inner crust are required, 

• 32 if the BPS (BPS) equation of state for outer crust and BBP equation of state for the 
inner crust are required, 

• 11 if the BPS (TMA) equation of state for the outer crust and polytropic equation of state 
for inner crust are required, and, 

• 12 if the BPS (BPS) equation of state for the outer crust and polytropic equation of state 
for inner crust are required. 

Equilibrium conditions 

• 1 for neutron star matter in general equilibrium, 

• 2 for neutron star matter in beta-equilibrium only, and, 

• 3 for pure neutron matter 

To illustrate the accuracy of the different subroutines for the different conventions the mass-radius 
relationship for the NL3 parameter set, assuming neutron star matter in general equilibrium and 
no crustal effects, is shown in Fig. IA.31 
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Figure A. 3: Mass-radius relationship for neutron star matter in general equilibrium, 
no crustal effects included, for the NL3 parameter set. 



APPENDIX B 
Solutions to integrals 



B.l Energy density 

Any expressions for the energy density of fermi-gas type model contains the integral 

1 f kp 



\ ! F dk k 2 y / k 2 + m 2 

7T 2 Jo 



where 

• k,F is the fermi momentum of the specific particle, and 

• m is the mass. 

This integral can be solved analytically by making the following substitution 

k 



> u, 

m 



and defining 

kp 
m 

thereby 



= Xp 



— - I dk k 2 V k 2 + m 2 = I duu 2 yu 2 + l 

71 Jo 71 Jo 

™ 4 f XF a 2 nr—, 

Jo 



Using integration by parts, the last integral reduces to: 

i 4 r x F _ , ™ 4 r 

* T JO 



% / du u Wu 2 + 1 = ^\ (24 + x F ) (1 + 4) 1/2 - sinh- 1 1 . (B.l) 

7T Jo O 71 " L 



B.2 Pressure 

As with the energy density, any expression for the pressure contains the integral 

1 f kF „ A; 4 
— 9 / are — . 

3?r Vfc 2 + m 2 
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By making the same substitution as with the energy density integral, the integral reduces to 



3tt 2 



dk 



Vk 2 + m? 



rri 

3^ 2 Jo 



kp/m 



du 



u 



m 



Vu 2 + 1 

4 



du ■ 



u 



37T 2 J y/y? + 1 



Which can be solved analytically as 



m 
3^2 



XF . u A m 4 1 



du 



(2x 3 F - 3x F ) (l + x 2 F ) 1/2 + 3sinh _1 (x F ) 



(B.2) 



B.3 MFT scalar meson field 

The expression <pQ ()9.20ap contains the integral 



1 

^ Jo 



kp , /c 2 m* 
dk 



Vk 2 + m* 



and as shown in Ref. [21] this integral can be expressed as 



•7T 2 / 

^ J 



A: 2 m* 



\/fc 2 + m* 2 vr 2 



2 , *2\ 1 / 2 *2 , 



kp I k F + m 



*2\V 2 ' 



. (B.3) 
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